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Abstract 

Use  of  Quantitative  Feedback  Theory  (QFT)  on  a  multiple-input  multiple- 
output  control  system  requires  certain  mathematical  properties  of  the  plant  ma¬ 
trix  of  system  transfer  functions  P(s).  In  general,  the  plant  matrix  P(s)  does  not 
possess  the  necessary  or  desired  mathematical  properties  for  the  QFT  design  to 
proceed.  A  frequency  sensitive  weighting  matrix  A(s)  is  used  to  transform  the 
plant  matrix  P{s)  into  an  equivalent  plant  matrix  Pe(s)  that  does  satisfy  QFT 
requirements.  In  matrix  notation,  the  relationship  between  the  equivalent  plant, 
plant,  and  weighting  matrices  is  Pe(s)  =  P(s) A(j)  .  This  thesis  identifies  the  nec¬ 
essary  and  desired  characteristics  of  the  equivalent  plant  Pt(s)  for  the  QFT  design 
process,  explains  the  use  of  the  Method  of  Specified  Outputs  which  generates  the 
frequency  sensitive  weighting  matrix  A(s),  and  calculates  several  A(s)  matrices 
for  a  3-input  2-output  lateral-directional  model  of  the  AFTI/F-16  aircraft.  For 
several  control  system  failures,  the  mathematical  structure  of  the  failed  equiva¬ 
lent  plant  matrices  Pe/(s)  is  examined  for  compliance  with  QFT  requirements.  A 
weighting  matrix  A(s)  is  found  that  produces  acceptable  equivalent  plant  matrices 
for  failures  of  down  to  0.01%  of  available  control  surface  deflections.  The  use  of 
the  software  packages  MATRIX*  and  MACSYMA  is  explained  as  applied  to  the 
Method  of  Specified  Outputs. 
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SELECTION  OF  A 

FREQUENCY  SENSITIVE  QFT  WEIGHTING  MATRIX  USING 
THE  METHOD  OF  SPECIFIED  OUTPUTS 

I.  Introduction 

Quantitative  Feedback  Theory  (QFT)  is  a  multiple-input  multiple-output 
(MIMO)  control  system  design  technique  that  accounts  for  uncertainties  inherent 
in  the  system’s  plant  state  space  model  [Hou87:l-lj.  Use  of  Quantitative  Feedback 
Theory  requires  certain  mathematical  properties  of  the  plant’s  state  space  model. 
The  system’s  plant  model  does  not  always  have  the  necessary  properties  for  the 
QFT  design  to  proceed.  The  goal  of  this  thesis  is  to  provide  the  tools  to  modify 
the  original  plant  matrix  into  a  form  that  can  be  used  in  the  QFT  design  process. 

For  a  MIMO  system,  the  plant’s  dynamics  can  be  described  as  a  matrix  of 
transfer  functions,  the  plant  matrix  P(s),  that  relate  the  m  control  inputs  UXS) 
to  the  l  system  outputs  T(s).  See  Figure  1.1.  Because  of  approximations  and 
measurement  limitations  in  determining  the  plant  model,  there  exists  an  enve¬ 
lope  of  uncertainty  for  the  plant’s  parameters  at  each  design  condition.  The  QFT 
technique,  which  requires  a  square  equivalent  plant  matrix  Pe(s),  attempts  to  de¬ 
termine  a  closed-loop  control  system  that  incorporates  a  cascade  compensator  G(s) 
and  an  input  filter  F(s)  that  provides  acceptable  system  performance  throughout 
the  range  of  plant  uncertainties.  See  Figure  1.2  for  the  QFT  compensated  MIMO 
system  [Hou87:3-7].  Note  that  in  Figure  1.2,  the  block  labeled  Pe(s)  is  the  square 
equivalent  plant  matrix  which  is  explained  below. 

As  stated  earlier,  the  plant’s  state  space  model  is  not  always  in  a  form  that 
allows  the  QFT  design  process  to  proceed.  QFT  design  requires  a  square  plant 
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Figure  1.1.  Uncompensated  Plant 


Figure  1.2.  QFT  MIMO  Control  Structure 
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Figure  1.3.  Equivalent  Plant 


matrix.  If  the  number  of  control  inputs  m  exceeds  the  number  of  system  outputs  /, 
then  an  “optimization”  problem  exists  for  exploiting  the  available  control  degrees 
of  freedom  [Hor88].  A  frequency  sensitive  weighting  matrix  A(a)  is  used  to  modify 
the  original  plant  matrix  in  order  to  exploit  the  available  degrees  of  control  freedom. 
Specifically,  the  system’s  plant  matrix  P{a)  is  postmultiplied  by  a  weighting  matrix 
A(s)  to  form  an  equivalent  plant  matrix  Pe(a).  See  Eq.  (1.1)  and  Figure  1.3.  Note 
that  Vi(s)  is  the  vector  of  desired  outputs. 

Pt(a)  =  P(s)  A(a)  (1.1) 

This  thesis  focuses  on  selection  of  a  frequency  sensitive  weighting  matrix 
A(s)  that  provides  the  necessary  mathematical  structure  to  the  plant  matrix  P(s) 
so  that  QFT  multiple-input  multiple-output  design  can  proceed.  The  Method  of 
Specified  Outputs  is  used  to  calculate  A(s).  Note  that  the  Method  of  Specified 
Outputs  technique  can  be  applied  to  other  multiple-input  multiple-output  control 
theories  that  use  a  weighting  matrix  to  transform  the  original  non-square  plant 
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matrix  into  a  square  equivalent  plant  matrix.  The  “squaring-down”  matrix  used 
by  Porter  is  an  example  [PM88:34-55]. 

1.1  Background 

Past  QFT  designs  have  used  a  weighting  matrix  A  of  constant,  real  elements 
S{j  in  order  to  transform  the  original  plant  matrix  P(s)  into  an  equivalent  plant 
matrix  Pe(s)  suitable  for  the  QFT  design  technique.  There  have  been  various 
methods  used  to  determine  the  magnitude  and  sign  for  each  A  element  6,j.  The 
first  attempts  in  determining  8{j  used  elements  of  0,  +1,  and  -1.  Engineering 
judgment  was  used  in  selecting  the  Sij’s.  One  technique  was  to  analyze  Bode  plots 
of  the  original  plant’s  transfer  functions  to  determine  which  control  inputs  had 
the  largest  influence  upon  the  system  outputs.  Later  QFT  designs  incorporated 
non-unity  Sij’s  but  they  were  still  real- valued  constants.  [Ham87:22] 

Last  year,  Lieutenants  Hamilton  and  Adams  began  using  frequency  depen¬ 
dent  elements  in  their  respective  digital  QFT  designs.  However,  LT  Hamil¬ 

ton’s  choice  of  the  structure  of  one  particular  (ij(s)  resulted  in  unacceptable  output 
and  he  was  forced  to  abandon  the  frequency  sensitive  ^;(s)  and  go  back  to  con¬ 
stant  valued  Sij’s.  [Ham87:39].  Also,  LT  Adams  optioned  to  use  constant  8tJ' s  in 
his  design.  [Ada88] 

1.2  The  Problem 

In  an  attempt  to  provide  QFT  designers  an  equivalent  plant  matrix  Pe(s) 
that  has  both  the  necessary  and  desired  structure  over  the  entire  range  of  original 
plant  P(s)  uncertainty,  a  frequency  sensitive  weighting  matrix  A(s)  is  sought. 
Additionally,  the  use  of  A(s)  on  the  original  plant  P(s)  is  actually  an  open-loop, 
cascade  compensation  of  the  original  open-loop  control  system  as  seen  in  Figure  1.1. 
Thus  the  use  of  a  weighting  matrix  can  improve  open-loop  performance  in  addition 
to  preparing  the  system  for  QFT  design.  The  primary  purpose  of  this  effort  is 
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to  find  a  A(s)  that  provides  an  equivalent  plant  Pe{s)  that  satisfies  the  QFT 
constraints. 


1.3  Scope 

To  achieve  a  solution  to  the  above  problem,  the  Method  of  Specified  Out¬ 
puts  is  used.  The  methodology  is  examined  and  a  lateral-directional  model  of  the 
AFTI/F-16  aircraft  is  used  as  an  example  of  the  method.  This  analysis  is  per¬ 
formed  only  in  the  continuous  Laplace  variable  s-domain.  Digital  considerations 
are  not  addressed. 

H  Assumptions 

It  is  assumed  that  the  QFT  designer  has  a  linear  state  space  model  of  the 
system  in  the  form 

x  =  Ax  +  Bu  (1.2) 

j/  =  Cx  4-  Du  (1.3) 

Where  x  =  state  variables,  u  =  control  inputs,  j/  =  system  outputs,  and  A,B,C, 
and  D  are  matrices  of  real,  constant  elements.  The  A,B,C ,  and  D  matrices  are 
only  valid  at  specific  design  conditions.  As  the  system’s  environment  changes,  so 
do  the  A,  J9,(7,and  D  matrices.  For  example,  if  the  system  modeled  is  an  aircraft, 
then  a  flight  condition  of  Mach  0.2  at  1,000  feet  altitude  produces  much  different 
A,  B,C,  and  D  matrices  than  a  flight  condition  of  Mach  1.5  at  30,000  feet  altitude. 

The  actual  system  being  modeled  determines  the  assumptions  needed  to 
transform  the  system’s  nonlinear,  coupled,  differential  equations  of  motion  into 
linear,  first  order  differential  equations  in  the  state  space  domain.  For  example, 
aircraft  state  space  model  assumptions  are  different  than  those  made  for  a  ship, 
a  robot,  or  a  missile.  For  the  Method  of  Specified  Outputs  of  Chapter  III,  it  is 
assumed  the  frequency  response  of  the  original  plant  matrix  P{s)  is  not  corrupted 
with  noise.  Thus,  the  coefficients  of  6ij{s)  is  uniquely  determined. 
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For  the  AFTI/F-16  aircraft  control  example  of  Chapter  V,  the  following 
assumptions  are  made: 

1.  The  aircraft  equations  of  motion  are  linearized  about  the  nominal  equilibrium 
(trim)  condition  of  Mach  0.9  at  20,000  feet. 

2.  A  small  perturbation  model  is  valid. 

3.  Input  commands  and  aircraft  response  do  not  violate  the  linear  model  as¬ 
sumptions. 

4.  The  mass  of  the  aircraft  remains  constant  during  the  commanded  maneuvers. 

5.  Aircraft  thrust  is  constant. 

6.  The  aircraft  is  a  rigid  body. 

7.  The  earth  is  an  inertial  reference  frame. 

8.  The  atmosphere  is  fixed  to  the  earth. 

9.  The  dimensional  stability  derivatives  for  flaperon  inputs  can  be  halved  to 
provide  derivatives  for  individual,  thus  independent,  surface  contributions. 

10.  The  dimensional  stability  derivatives  are  constant  for  the  given  flight  condi¬ 
tion. 

11.  Control  system  failures  can  be  described  by  scaling  columns  of  the  control 
matrix  B  of  the  state  space  model  Eq.  (1.2).  Scaling  refers  to  a  positive 
scalar  less  than  one  multiplying  an  appropriate  column  of  the  B  matrix.  This 
assumption  means  the  aircraft’s  stability  derivatives  taken  with  respect  to  the 
states,  which  are  contained  in  the  A  and  C  matrices  of  Eqs.  (1.2)  and  (1.3), 
do  not  appreciably  change  when  a  control  surface  is  failed  as  either  frozen  in 
an  arbitrary  position  or  is  unable  to  achieve  maximum  deflection.  With  this 
assumption,  many  control  system  failures  can  be  modeled. 

12.  Laplace  transforms  exists  for  all  signals  and  all  system  elements. 

See  [Bar83j  and  [Sch86j. 


1-6 


1.5  Materials  and  Equipment 

The  analysis  is  performed  using  MATRIX*  software  on  a  VAX  11-780  digital 
computer.  MATRIX*  is  a  programmable  matrix  calculator  with  graphics  and  it 
has  extensive  control  system  analysis  capabilities.  [MAT86] 

Also,  MACSYMA  software  on  a  VAX  11-785  digital  computer  is  used  to 
perform  symbolic  (as  opposed  to  numerical)  matrix  multiplication,  inversion,  and 
determinants.  [VAX85] 

1.6  Presentation 

Chapter  II  describes  the  QFT  constraints  on  the  equivalent  plant  matrix. 
These  in  turn  constrain  the  degrees  of  freedom  for  choosing  the  weighting  matrix 
A (s)  that  transforms  the  plant  matrix  P(s)  into  the  equivalent  plant  matrix  Pe(s). 

Chapter  III  describes  the  Method  of  Specified  Outputs.  This  method  is  used 
to  determine  the  frequency  sensitive  weighting  matrix  A(s). 

Chapter  IV  describes  the  selection  criteria  for  the  QFT  weighting  matrix 
based  upon  the  QFT  constraints  of  Chapter  II  and  the  Method  of  Specified  Outputs 
of  Chapter  III. 

Chapter  V  explores  an  example  of  the  Method  of  Specified  Outputs  tech¬ 
nique  using  a  lateral-directional  model  of  the  AFTI/F-16.  Analysis  and  results  are 
presented. 

Chapter  VI  summarizes  conclusions  of  this  work  and  recommendations  for 
further  investigation. 

Appendix  A  provides  examples  of  MATRIX*  command  files  useful  for  the 
numerical  aspect  of  the  Method  of  Specified  Outputs. 

Appendix  B  provides  examples  of  MACSYMA  commands  useful  for  the  sym¬ 
bolic  matrix  manipulation  aspect  of  the  Method  of  Specified  Outputs. 
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II.  QFT  and  the  Equivalent  Plant  Matrix 

2.1  Introduction 

Use  of  Quantitative  Feedback  Theory  requires  the  plant  matrix  of  transfer 
functions  to  have  certain  mathematical  properties.  In  general,  the  plant  does  not 
possess  the  necessary  and  desired  properties  for  the  QFT  design  technique  to  be 
applied.  Therefore,  the  plant’s  state  space  model  is  modified  into  the  desired  form 
by  open-loop  weighting  matrix  compensation.  See  Figure  1.3. 

In  the  past,  researchers  used  a  single  weighting  matrix  of  constant,  real  el¬ 
ements  to  perform  the  necessary  mathematical  modification  of  the  plant  matrix 
[Ham87:22].  They  had  some  success.  However,  this  thesis  is  an  attempt  to  investi¬ 
gate  a  technique  to  find  a  single,  frequency  sensitive  weighting  matrix  A(s)  which 
provides  an  equivalent  plant  transfer  function  matrix  Pe(s)  with  the  necessary  and 
desired  mathematical  properties  over  the  entire  range  of  plant  uncertainties. 

2.2  Plant  Matrix  P(s) 

The  plant  matrix  P(s)  is  the  transfer  function  matrix  that  relates  all  system 
inputs  to  all  system  outputs.  The  system  is  modeled  as  a  set  of  linear,  first  order, 
vector  differential  equations 

x  =  Ax  4-  Bu  (2.1 ) 

y  =  Cx  +  Du  (2.2) 

where  x  is  the  vector  of  state  variables,  u  is  the  vector  of  control  inputs,  y  is  the 
vector  of  system  outputs,  and  the  A,  B,C ,  and  D  matrices  contain  real  constants. 
The  elements  of  the  state,  input,  and  output  vectors  are  functions  of  time.  Let  m  be 
the  number  of  control  inputs  and  l  be  the  number  of  system  outputs.  The  vector 
u  has  dimension  mxl  and  the  vector  y  has  dimension  fxl.  Taking  the  Laplace 
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transform  of  these  equations  and  assuming  zero  initial  conditions  gives 


£_={sl  -  A)~lBU_  (2.3) 

r=C£+Z?If  (2.4) 

where  £,  H,  and  Y_  are  the  Laplace  transforms  of  x(t),  u(t)  and  t/(t)  respectively. 
Solving  for  Y(s)  in  terms  of  A,B,C,D  and  U(s)  gives 

11=  [C{sl~  A)~1B  +  d]U_  (2.5) 

Let  the  plant  matrix  P(s)  be  defined  as 

P{s)  3  C{sl-  A)~XB  +  D  (2.6) 

Note  that  P(s)  is  an  uncompensated,  open-loop  control  system  with  m  inputs  and  l 
outputs  as  seen  in  Figure  1.1.  The  plant  P(s)  is  a  ixm  matrix  of  transfer  functions. 

The  plant  matrix  P(s)  embodies  the  multiple-input  multiple-output  (MIMO) 
system  model  as  seen  in  Figure  1.1  and  in  Eq.  (2.5).  A  different  plant  matrix  P(s) 
exists  for  each  different  environment  or  system  failure  encountered.  For  example, 
if  the  system  being  modeled  is  an  aircraft,  a  different  P(s)  exists  for  each  flight 
condition  and  for  each  control  system  failure  [i.e.  for  any  possible  combination  of 
the  plant  parameters  contained  in  A,  B,C,  or  D  of  Eq.  (2.6)]. 

For  the  entire  range  of  P{s)' s  at  a  given  design  condition,  there  exists  math¬ 
ematical  properties  that  are  either  necessary  or  very  desirable  for  the  MIMO  QFT 
design  technique.  In  general,  P{s)  does  not  possess  the  necessary  or  desired  prop¬ 
erties.  To  solve  this  problem,  a  weighting  matrix  A(s)  is  used  to  transform  P[s)  in 
to  an  equivalent  plant  matrix  Pe(s)  that  does  satisfy  the  QFT  requirements.  It  is 
hoped  that  all  P(s)'s  within  the  region  of  plant  uncertainty  can  be  made  to  satisfy 
the  QFT  requirements  with  the  use  of  a  single  A(s). 
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2.3  Equivalent  Plant  Pe(s)  and  Weighting  Matrix  A(s) 

The  equivalent  plant  is  obtained  by  post-multiplying  the  original  nominal 
plant  matrix  by  a  weighting  matrix  A(s). 

Pe(s)  =  P(s)A(3)  (2.7) 

The  basis  for  Eq.  (2.7)  can  be  seen  from  the  following  observations  of  Figure  1.3: 
The  original,  uncompensated,  input-output  relationship  is  given  by 

Z  =  P(s)LL  (2.8) 

By  introducing  a  weighting  matrix  A(s)  cascaded  with  the  original  plant  matrix 
as  in  Figure  1.3,  the  input-output  relation  for  the  A(s)  matrix  is 

[L  =  A(s)V  (2.9) 

where  the  vector  V_  is  the  Laplace  transform  of  the  desired  output  vector  v(t). 
Substituting  Eq.  (2.9)  into  Eq.  (2.8)  gives  the  equivalent  plant  input-output  rela¬ 
tionship  of  Figure  1.3 

r  =  P(s)A(s)F  (2.10) 

The  equivalent  plant  Pe(s)  is  defined  as 

P<(s)  =  P(*)  A(*)  (2.11) 

Thus  the  overall  input-output  relationship  Eq.  (2.10)  can  be  written 

Y.  =  Pt(*)V  (2-12) 

Recall  that  V_  is  the  vector  of  Laplace  transformed  desired  outputs  and  Y  is 
the  vector  of  Laplace  transformed  actual  outputs.  Therefore,  defining  P'(s)  such 
that  Pt(s)  =  P(s)A(s)  provides  the  designer  a  mechanism  to  modify  the  system 
plant  model  into  a  form  useful  for  QFT  design.  The  A(s)  matrix  is  the  mechanism 
to  achieve  the  useful  form. 
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Note  that  the  use  of  A(s)  is  essentially  an  open-loop  matrix  compensation  of 
the  original  MIMO  system.  Improvement  in  system  performance  is  expected  with 
proper  choice  of  the  weighting  matrix  A(s)  in  addition  to  preparing  the  system 
model  for  use  by  the  Quantitative  Feedback  Theory’s  technique  [Lew88].  Of  course, 
the  number  one  priority  of  the  weighting  matrix  A(a)  is  for  QFT  purposes. 

2. 4  Characteristics  of  the  Equivalent  Plant  Pe(a) 

i.f.l  QFT  Requirements  on  the  Equivalent  Plant  This  section  deals  with 
the  required  mathematical  structure  of  the  equivalent  plant  matrix  Pt(s).  The 
justification  for  these  requirements  come  from  the  Multiple-Input  Single-Output 
(MISO)  Equivalent  Method  used  in  MIMO  QFT  designs  [Hou87:3-9].  Multiple- 
input  multiple-output  QFT  designers  reduce  their  MIMO  systems  to  a  set  of  equiv¬ 
alent  MISO  systems.  To  use  the  MISO  equivalent  models,  the  following  conditions 
must  be  met. 

Necessary  Condition  1  The  equivalent  plant  matrix  Pr(s)  must  be  non-singular 
for  both  the  nominal  plant  and  for  the  range  of  plant  uncertainties  [Hou87:3-15]. 
Thus  P,(s)  must  be  controllable  [Hor88j.  The  equivalent  plant  matrix  must  be 
square  and  the  determinant  of  Pe(s)  must  be  non-zero  for  any  possible  combina¬ 
tions  of  plant  parameters  [Hou87:3-15j. 

A  square  equivalent  plant  matrix  requirement  means  that  the  number  of  con¬ 
trol  inputs  equals  the  number  of  system  outputs.  In  any  multiple-input  multiple- 
output  system  with  m  inputs,  there  are  at  most  rn  outputs  which  can  be  inde¬ 
pendently  controlled.  [Hou87:3-2].  Thus  the  dimensions  of  Pe(s)  will  be  Ixl  where 
l  <  rn. 

It  is  easy  to  find  an  example  where  the  original  plant  P(s)  is  not  square. 
Consider  a  control  problem  where  there  are  more  control  inputs  available  than 
outputs  that  the  designer  wishes  to  control  (i.e.  P(s)  has  more  columns  than 
rows).  For  example,  if  the  system  being  modeled  is  an  aircraft,  there  are  often 
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more  control  surfaces  than  aircraft  output  variables  that  the  designer  is  interested 
in  controlling.  For  example,  in  the  lateral- directional  model  of  an  aircraft,  one 
may  desire  to  control  sideslip  and  roll  rate  but  the  available  independent  control 
surfaces  are  left  and  right  flaperons,  left  and  right  differential  horizontal  tail,  and 
rudder.  In  this  case,  there  are  5  control  inputs  and  2  system  outputs  and  the  plant 
matrix  P(s)  is  2x5.  In  order  for  Pe(s)  to  be  2x2,  then  A (s)  must  be  5x2. 

Necessary  Condition  2,  Only  for  the  2x2  Case  For  an  equivalent  plant  ma¬ 
trix  with  2  inputs  and  2  outputs 


PM 


P*ll  Pe  n 
Pe  21  Pen 


(2.13) 


the  diagonal  dominance  condition  for  the  2x2  equivalent  plant  states  that  as  u>  — +  oo 
then 


P<  u  Pen 


> 


P*ijP*ji 


(2.14) 


See  [Hou87:3-I7 


Necessary  Condition  3,  Only  for  the  3x3  Case  The  effective  plant  matrix  for 
the  3x3  case  is 


Pen 

Pen 

Pen 

P  eji 

Pe  n 

Pe  n 

Pe„ 

Pe  jj 

Pe„ 

(2.15) 


and  its  inverse  is 


P(s)-1 


Pen 

Pen 

Pen 

p;„ 

Pen 

P^n 

Pe'„ 

P«'n 

P' 53 

(2.16) 


If  the  QFT  designer  is  using  the  Single-Loop  Equivalents  technique  of  [Hou87:3-19] 
on  a  3x3  plant,  then  the  following  condition  must  hold:  as  w  ->  oc 


Pe  1 1  P«jj  Pe,,  —  Pen  ( ]Pe,i  Pe,,  Pe„Pe,i  |)  |p*n  (  Pen  Pe,i  Pe3i  Pe„  )  ~  Pei  i  P  e2,  Pe, 


(2.17) 
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If  the  QFT  designer  is  using  the  improved  method  instead  of  the  Single- Loop 
Equivalents  technique,  then  the  following  condition  must  hold:  as  u>  — ►  oo  then 
det  [|rj  must  not  change  sign  where  P0  is  the  nominal  plant  [Hor88j.  See  [Hor79] 
for  higher  order  plants. 

2-4-2  Desired  Characteristics  of  the  System  Model  This  section  deals  with 
the  desired  characteristics  of  the  equivalent  plant  matrix  Pt{s). 

Desired  Characteristic  1  The  determinant  of  the  equivalent  plant  matrix  Pe(s) 
should  preferably  be  minimum  phase  for  the  nominal  equivalent  plant,  and  for  all 
possible  combinations  of  control  system  failures.  Thus  it  is  desired  that  det[Pe(s)] 
be  minimum  phase  for  the  entire  range  of  plant  uncertainties.  That  is,  since  the 
elements  of  Pt(s)  are  quotients  of  polynomials  in  the  Laplace  variable  s,  then  the 
zeros  of  the  determinant  of  P^s),  det[Pe(j)j,  should  he  in  the  left-half  s-plane  for 
both  the  healthy  and  failed  plants.  If  zeros  of  det[Pe(a)]  do  he  in  the  right-half 
s-plane,  then  they  should  preferably  be  far  into  the  r.h.p.  so  as  to  minimize  their 
effect  [Hor88]. 

A  sufficient  condition  for  det[P«(s)]  to  be  made  minimum  phase  is  that  at 
least  one  of  the  Ixl  minors  of  P(  s)  be  minimum  phase  over  the  entire  range  of  plant 
parameter  variations  [Hou87:5-19]  and  [Sch86:ll].  This  conclusion  comes  from  the 
Binet-C’auchy  formula  [Lan69]. 

As  noted  previously,  this  thesis  does  not  consider  the  digital  QFT  design 
technique.  The  digital  QFT  designer  should  realize  that  sampled  systems  are 
inherently  non-minimum  phase  [Hou87:E-8]. 

Desired  Characteristic  2  The  equivalent  plant  matrix  should  be  as  diagonal 
as  possible  [Hor88].  If  Pe{s)  is  diagonal,  then  no  cross-coupling  occurs  between 
control  inputs  and  system  outputs.  That  is,  one  input  effects  only  one  output. 
Generally,  however,  all  failed  equivalent  plants  Ptf($)  are  not  diagonal. 

Desired  Characteristic  3  For  a  stable  original  plant  P(s),  the  response  of  the 
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2-6 


equivalent  plant  Pe(s)  should  be  reasonable  and  stable  for  reasonable  control  inputs 
for  the  entire  range  of  plant  parameter  variations.  For  stability,  the  poles  of  the 
equivalent  plant  matrix  (i.e.  roots  of  the  denominator  polynomial  for  each  plant 
element)  should  lie  in  the  left-half  j-plane.  A  reasonable  response  means  that  the 
outputs  are  not  excessive  provided  that  the  inputs  are  not  unusually  large  and 
that  both  the  input  and  output  limits  are  not  exceeded.  The  reasonableness  of  a 
response  is  a  judgement  made  by  the  QFT  designer.  The  designer  must  consider 
the  maximum  and  minimum  input  control  limits  and  output  response  limits  of 
the  system.  For  example,  if  the  system  being  modeled  is  an  aircraft  and  if  a  step 
rudder  input  of  5  degrees  causes  the  sideslip  angle  to  be  greater  than  95  degrees, 
then  that  scenario  would  fail  the  reasonability  test.  For  the  same  aircraft,  if  an 
aileron  deflection  of  200  degrees  were  required  to  cause  a  roll  rate  of  1  degree  per 
second,  that  scenario  would  also  fail  the  reasonability  test. 

A  reasonable  response  implies  that  the  elements  of  Pt  must  be  realizable. 
That  is,  each  transfer  function  of  Pe(s)  must  have  at  least  as  many  poles  as  zeros. 

The  remainder  of  this  thesis  focuses  on  the  selection  of  a  frequency  sensitive 
weighting  matrix  A (s)  by  use  of  the  Method  of  Specified  Outputs.  For  the  example 
of  Chapter  5,  the  resulting  equivalent  plant  Pt(s)  is  in  the  correct  form  for  the  QFT 
design  to  proceed. 
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>  ///.  T7te  Method  of  Specified  Outputs 

3.1  Introduction 

The  Method  of  Specified  Outputs  allows  calculation  of  a  frequency  sensitive 
I 

weighting  matrix  A(s)  for  use  with  the  Quantitative  Feedback  Theory  design  pro¬ 
cedure.  The  Method  of  Specified  Outputs  was  developed  by  Lieutenant  Colonel 
Zdzilaw  H.  Lewantowicz  [Ham87j.  As  noted  in  Chapter  II,  the  method  can  be 
*  thought  of  as  an  open-loop  compensation  of  the  system.  The  goal  is  to  find  a 

single  weighting  matrix  A(s)  which  provides  an  acceptable  Pe{s)  for  not  only  the 
healthy  P(s)  but  also  for  the  entire  range  of  failed  plant  matrices  P/(s). 
l 

3.2  The  Method  of  Specified  Outputs 

The  Method  of  Specified  Outputs  is  based  upon  Eq.  (3.1) 

1  P,W  =  P(a)A(a)  (.5.1) 

The  key  to  the  method  is  specifying  the  equivalent  plant  matrix  Pe(s)  that  satisfies 

the  QFT  constraints  described  in  Chapter  II.  If  Pe(s)  and  P{s)  are  given,  then 
I 

A (s)  can  be  determined  from  linear  algebra  transformations  discussed  below. 

CASE  1:  If  the  original  plant  has  more  control  inputs  than  system  outputs, 
then  P{s)  has  more  columns  than  rows  and  A(s)  can  be  defined  as  the  minimum 

> 

norm  solution  to  Eq.  (3.1)  provided  that  P{s)  has  linearly  independent  rows.  The 
minimum  norm  solution  to  Eq.  (3.1)  is 

A(3)  =  P(s)rlP(a)P(i)T]-1P.(S)  (3.2) 

CASE  2:  If  the  original  plant  P{s)  has  the  same  number  of  inputs  as  outputs 
and  if  P{s)~l  exists  but  P{s)  does  not  meet  the  QFT  constraints  of  Chapter  II, 
then  A  is  given  as 
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A(,)  =  P(s)-lP'(s) 


(3-3) 


CASE  3:  If  the  original  plant  P(a)  has  fewer  control  inputs  than  system 
outputs,  some  of  the  outputs  will  be  uncontrollable.  In  this  case,  an  exact  solution 
for  A(s)  is  not  possible  but  a  least  squares  solution  is  given  by  Eq.  (3.4)  provided 
that  P(s)  has  linearly  independent  columns. 

A(»)  =  [■P(,)rP(.)]-,J>(,)T*W  (3.4) 

3.2.1  The  Transformation  of  P(s)  into  Pg(a)  The  A  (a)  matrix  transforms 
the  P(.j)  matrix  into  the  P€(s)  matrix.  The  columns  of  the  A(s)  matrix  “steer” 
linear  combinations  of  the  columns  of  the  P{s)  matrix  into  the  columns  of  the 
Pe(a)  matrix  [Lew88].  For  the  minimum  norm  and  exact  solutions  of  cases  1  and  2 
above,  the  transformation  of  P(s)  into  P«(a)  space  is  exact.  For  the  least  squares 
solution  of  case  3  above,  the  P  matrix  can  not  be  transformed  into  Pt  space  but 
the  error  of  the  projection  of  P  into  Pt  space  can  be  minimized.  The  following 
section  describes  Eq.  (3.2)  thru  Eq.  (3.4).  The  selection  of  Pt(s)  is  discussed  in 
Chapter  IV. 

3. 2. 1.1  The  Minimum  Norm  Solution  For  the  case  where  P(s)  has 
more  columns  than  rows,  the  minimum  norm  solution  for  A(s)  is  given  by  Eq.  (3.5). 
When  P(s)  has  more  columns  than  rows,  this  is  equivalent  to  saying  that  there  are 
more  control  inputs  u  than  system  outputs  y.  This  is  an  underconstrained  matrix 
algebra  problem  and  there  are  infinitely  many  exact  solutions  for  A(s),  provided 
that  the  rank  of  P(s)P(s)T  is  as  large  as  possible.  The  minimum  norm  solution 
for  A(s)  is 

A(s)  =  P(s)r(P(s)P(s)rj-1Pe(s)  (3.5) 
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The  expression  P(s)T(P(s)P(s)T)~l  is  called  the  right  inverse  of  P(a)  and 
is  thus  a  pseudoinverse  of  P(s)  [Gel74:19].  If  P(s)  is  post-multiplied  by  its  right 
inverse,  then  the  identity  matrix  results.  If  the  minimum  norm  solution  for  A  is 
substituted  into  Eq.  (3.1),  then  the  following  identity  results 

PM  =  PU)P(’f{P(‘)P(‘)T\-'PM  (3.6) 


or 


Pe(s)  =  Pe(s) 


(3.7) 


The  necessary  condition  for  the  minimum  norm  solution  for  A(s)  to  exist  is 


if 


rank{PPT)  =  l 


(3.8) 


where  /  =  row  dimension  of  P. 


Eq.  (3.8)  says  P(s)  must  have  linearly  independent  rows  and  that  rank{P PT) 
is  as  large  as  possible.  If  Eq.  (3.8)  holds,  then  a  unique  [PPT)~l  exists,  the  right 
inverse  of  P{s)  exists,  and  the  minimum  norm  solution  for  A(s)  exists  [Str88:96]. 
Note  that  ( PPT )  is  an  Ixl  square  symmetric  matrix. 

The  minimum  norm  solution  is  only  one  of  infinitely  many  possible  solutions 
for  the  underconstrained  case  where  P(s)  has  more  columns  than  rows.  The  term 
“minimum  norm”  comes  from  the  fact  that  the  right  inverse  (i.e.  PT[P Pr]_1) 
minimizes  the  norm  of  each  column  of  A(s).  That  is,  each  column  of  A(s)  is  a 
solution  vector  of  minimum  length  [Gel74:19]. 

This  minimum  norm  solution  for  A(s)  is  used  in  the  example  of  Chapter  V. 


3. 2. 1.2  The  P~x  Solution  For  the  case  where  P(s)  has  the  same  num¬ 
ber  of  columns  and  rows,  and  if  rank(P)  =  l  =  m,  then  a  unique  P(s)~1  exists. 
Only  one  solution  for  A (s)  exists  and  is  given  by 


A(s)  =  P(s)-1Pf(s) 


(3.9) 
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The  QFT  designer  can  use  Eq.  (3.9)  if  P(s)  does  not  possess  the  desired 
plant  matrix  properties  described  in  Chapter  II.  By  properly  specifying  Pe  and 
using  Eq.  (3.9),  the  desired  equivalent  plant  matrix  properties  can  be  achieved  for 
at  least  the  healthy  plant. 

3.Z.1.3  Least  Squares  Solution  For  the  case  where  P(s)  has  more  rows 
than  columns  and  P(s)  has  linearly  independent  columns,  then  the  least  squares 
solution  for  A (s)  is 

A(»)  =  lPMTP(,)r‘P(s)TP,(s)  (3.10) 

For  Eq.  (3.10)  to  have  a  valid  solution,  then  the  following  condition  must  be 

met 

rank(PT  P)  =  m  (3-11) 

where  m  =  column  dimension  of  P.  Eq.  (3.11)  says  P  must  have  linearly  indepen¬ 
dent  columns  and  that  ranlc(PT P)  is  as  large  as  possible  [Str88:15|. 

When  P(s)  has  more  rows  than  columns,  this  is  equivalent  to  saying  there  are 
more  system  outputs  t/  than  control  inputs  u.  Obviously,  the  excess  of  outputs  t/ 
over  the  inputs  u  can  never  be  controlled.  This  is  an  overconstrained  matrix  algebra 
problem  and  no  exact  solution  exits.  Eq.  (3.10)  guarantees  the  error  between  Pt(s) 
and  P(s)A(s)  is  minimized  in  the  least  squared  sense. 

Note  that  PT P  is  mxm  square  and  symmetric.  If  the  condition  of  Eq.  (3.11) 
holds,  then  (PTP)~l  exists,  the  left  inverse  of  P  exists,  and  Eq.  (3.10)  applies. 
The  left  inverse  of  P  is  defined  as  (PT P)_1  PT  and  is  a  pseudoinverse  of  P.  If  P  is 
premultiplied  by  its  left  inverse,  the  identity  matrix  results,  i.e., 

( PTP)~lPTP  =  I  (3.12) 

Substituting  Eq.  (3.10)  into  Eq.  (3.1)  gives 

Pt,opprx  =  P(PTP)-lPTPt  (3.13) 
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Note  the  left  hand  side  of  Eq.  (3.13)  is  denoted  Pe,apprx  because  it  is  not 
exactly  the  specified  Pe  since  no  exact  solution  exists.  The  matrix  P{PT P)~x PT 
projects  Pe  onto  the  column  space  of  P  and  Pt,apprx  is  that  projection  [Str88:158]. 

3.2.2  Numerical  Determination  of  A  (a)  The  above  section  implies  that  the 
calculation  of  A(s)  is  performed  symbolically.  That  is,  for  the  minimum  norm  so¬ 
lution  A  =  PT(PPT)~1Pe  ,  there  is  a  matrix  product  calculation  PPT ,  a  matrix 
inversion  ( PPT )  l,  and  two  more  matrix  products  PT(PPT)  1  and  PT(PPT)  xPe 
before  A(s)  is  determined.  Recalling  that  each  element  of  P( s)  is  a  transfer  func¬ 
tion  in  3  (i.e.  the  quotient  of  two  polynomials  in  the  Laplace  variable  3)  it  is  clear 
that  the  symbolic  calculation  of  A (3)  is  time  consuming  if  performed  by  hand 
calculation.  It  is  difficult  for  even  a  1x2  plant  with  a  denominator  polynomial  of 
degree  3  and  the  numerator  polynomials  of  degree  2.  To  speed  up  the  calcula¬ 
tion  of  A(j),  the  QFT  designer  can  use  a  symbolic  manipulator  like  MACSYMA 
[VAX85]  or  the  QFT  designer  can  determine  A(a)  numerically  at  selected  complex 
frequencies  3*  =  j u>fc  and  determine  the  denominator  and  numerator  polynomial 
coefficients  using  a  systems  identification  scheme  [Ham87].  Note  that  numerically, 
A(jwfc)  is  calculated  exactly  because  both  P{jw\,)  and  Fe(ju;*.)  are  known  exactly 
by  evaluating  F(s)  and  F,(s)  at  3  =  jwj,. 

Both  conceptual  approaches  for  determining  A(s)  are  examined  for  the  exam¬ 
ple  of  Chapter  V.  It  is  not  possible  to  sufficiently  simplify  the  symbolic  MACSYMA 
expression  for  A(s)  into  a  manageable  form.  Therefore,  the  numerical  approach 
is  chosen  for  the  example  of  Chapter  V.  Once  the  transfer  functions  of  A(s)  are 
determined  numerically  using  MATRIX*,  then  MACSYMA  is  used  to  perform  the 
matrix  product  Pe(s )  =  P(s)A(s)  . 

Determining  A(s)  numerically  is  essentially  a  mathematical  systems  identi¬ 
fication  problem.  Two  such  system  identification  techniques  are  considered.  One 
technique  was  devised  by  Lieutenant  Colonel  Z.  Lewantowicz  and  the  other  by  Lieu- 


tenant  Colonel  W.  Baker  both  of  the  Air  Force  Institute  of  Technology.  The  scheme 
devised  by  LTC  Baker  seems  well  suited  for  this  MIMO  application  but  there  was 
insufficient  time  to  investigate  this  technique  thoroughly.  See  reference  [Bak80]  for 
a  MIMO  systems  identification  technique  and  references  [Lev59]  and  [PayTOj  for  a 
SISO  systems  identification  technique.  The  remainder  of  this  section  discusses  the 
determination  of  A (s)  using  LTC  Lewantowicz’s  systems  identification  technique. 

Evaluate  the  appropriate  A(s)  solution  of  Section  3.2.1  at  a  number  of  com¬ 
plex  frequencies  s  =  j u>.  The  minimum  number  of  necessary  frequency  evaluations 
/  is  discussed  below.  Consider  an  individual  element  of  the  numerically  computed 
weighting  ipatrix  6ij(jcj).  Each  element  is  a  complex  number  which  is  the  evalua¬ 
tion  of  a  transfer  function  at  a  specific  frequency.  The  element’s  transfer  function 
is  of  the  form 

'AJU)  br(ju,y  +  +  . . .  +  bi(juj)  +  bo  1  } 

where  q  is  the  degree  of  the  numerator  polynomial  and  r  is  the  degree  of  the  denom¬ 
inator  polynomial.  The  goal  is  to  determine  the  real-valued  constant  coefficients 
a,,  a,_ j,  •  •  • ,  oj,  ao,  br,  6r-i>  •  •  • ,  &i,  60  for  each  element  Sij  of  the  weighting  matrix 
A (s).  One  reason  for  calculating  these  coefficients  is  the  desire  to  determine  the 
time  histories  of  the  outputs  of  A( s).  For  the  transfer  function  to  be  realizable  (i.e. 
proper),  then  q  <  r.  Note  that  the  outputs  of  A(j)  are  the  inputs  to  the  original 
plant  P(s).  See  Eq.  (2.9).  The  method  used  to  calculate  A(^)  imposes  no  path 
constraints  on  the  time  histories  of  the  outputs  of  A  and  thus  should  be  checked 
for  reasonableness. 

The  following  technique  determines  the  coefficients  for  each  ^(s).  Once  all 
£jj(s)’s  have  been  determined,  a  common  denominator  polynomial  for  the  A(s) 
matrix  can  be  found.  Again  note  that  when  Sij(ju)  is  evaluated  at  a  particular 
complex  frequency  ju)k,  then  Sij(ju>it)  is  a  complex  number.  In  matrix  notation, 
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the  scalar  Eq.  (3.14)  becomes 


(3.15) 

Further  manipulation  of  Eq.  (3.15)  yields  a  matrix  equation  that  evaluates 
to  the  scalar  zero 


(3.16) 
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Define  c,j  as  the  vector  of  coefficients  that  define  the  transfer  function  Sij(juj) 


£.i  = 


(3.17) 


Cq—i 

CLi 

Co 

br 

br- 1 

bi 
bo 

Note  that  dj  has  the  dimension  (g  +  r  +  2]xl. 

The  objective  of  the  system  identification  technique  is  to  identify  each  ele¬ 
ment  of  This  can  be  done  because  can  be  calculated  exactly  at  any 

frequency  jui j,.  To  use  LTC  Lewantowicz’s  technique,  evaluate  Eq.  (3.16)  at  / 


uencies 

where  /  >  (q  + 

r +  2) 

(M)’ 

(j^i)*"1  ••• 

(M)  1 

-Hjvi) 

(j“i  )r 

...  i 

(-7^2  )" 

••• 

(Jvt)  1 

—  b(j  U>2) 

(jv2y 

...  i 

(M)q 

(jW/),_1  ••• 

(j^f)  1 

-Hjvf) 

{jvfY 

...  i 

a,-i 

a, 

Co 

br 
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Note  that  Eq.  (3.18)  has  the  form  Fc.j  =  Q.  Let  F  be  the  rectangular 
}x[q  +  r  -f-  2\  matrix  of  complex  numbers  that  premultiplies  cx}.  1  herefore,  ctJ  lies 
in  the  nullspace  of  F  where  F  must  be  singular  for  a  nullspace  to  exist.  Note  that 


if  F  is  non-singular,  then  the  QFT  designer  should  evaluate  the  numerical  solution 
to  A(s)  at  more  distinct  frequencies  so  that  F  is  singular  and  a  nulispace  exists. 

The  vector  defines  only  one  transfer  function  of  the  weighting  matrix  A(s). 
Since  A(a)  has  dimension  mxi,  then  the  procedure  to  find  c^  must  be  repeated 
m  x  l  times.  There  are  two  ways  of  determining  the  first  is  to  transform  F 
into  modified  Hermite  normal  form.  The  second  way  is  to  perform  a  singular  value 
decomposition  on  F. 

Modified  Hermite  Normal  Form  Technique.  Let  F  have  at  least  as  many  rows 
as  columns.  Let  F  be  the  modified  Hermite  normal  form  of  F.  The  modified 
Hermite  normal  form  of  a  matrix  is  upper  triangular  (i.e.  all  elements  below 
the  main  diagonal  are  zero),  has  either  ones  or  zeros  on  the  main  diagonal,  and 
any  column  with  unity  on  the  diagonal  is  otherwise  zero  [Rei83:330-331].  Note 
that  a  matrix  in  modified  Hermite  normal  form  is  sometimes  referred  to  as  row 
reduced  echelon  form  [Str88:77|.  Replace  any  zero  on  the  main  diagonal  of  F 
with  -1.  Let  Z  be  a  matrix  composed  of  the  first  q  +  r  +  2  elements  of  any 
column  of  F  that  has  a  —1  on  the  main  diagonal.  The  columns  of  Z  span  the 
nulispace  of  F  and  therefore  any  column  of  Z  could  be  chosen  to  be  the  vector 
of  coefficients  c .  However,  care  should  be  taken  when  selecting  which  column 
of  Z  should  be  chosen  as  c,j.  F  is  a  complex  matrix  and,  in  general,  a  set  of 
complex  vectors  lie  in  the  nulispace  of  F.  However,  in  order  for  the  transfer 
function  8,j  to  have  physical  meaning,  the  elements  of  cxj  must  be  real.  Therefore, 
chose  the  column  of  Z  which  has  the  negligible  imaginary  components  and  use  the 
real  components  as  c^.  To  check  that  the  real-valued  does  indeed  he  in  the 
nulispace  of  F,  evaluate  the  matrix  product  Fc,;.  This  product  should  be  an  /xl 
vector  of  elements  with  each  element’s  magnitude  less  than  about  10-10.  Another 
useful  check  is  to  numerically  evaluate  the  transfer  function  defined  by  ctJ  and 
compare  those  complex  numbers  with  the  appropriate  elements  of  the  numerically 
determined  A(s)  solution  of  Section  3.2.1. 
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MATRIX*  lias  a  row  reduced  echelon  routine  called  RREF(matrix)  [MAT86]. 
See  Appendix  A. 

Singular  Value  Decomposition  Technique.  Any  fx[q  +  r  -f  2]  matrix  F  can  be 
factored  into 

F  =  Qli:Q2  (3.19) 

where  the  £  singular  values  of  F  fill  the  first  £  places  on  the  main  diagonal  of 
the  fx[q  +  r  +  2]  dimensioned  E  matrix,  the  columns  of  the  fxf  dimensioned  Q1 
matrix  are  the  eigenvectors  of  FFT ,  and  the  columns  of  the  [g  -f  r  4-  2]x[g  +  r  -f  2] 
dimensioned  Q2  matrix  are  the  eigenvectors  of  FTF.  The  last  [g  +  r-f  2]  —  £  columns 
of  Qi  he  in  the  nullspace  of  F  (Str88:443j. 

Again,  any  real  vector  that  lies  in  the  nullspace  of  F  can  be  used  as  the  vector 
of  coefficients  of  the  transfer  function  Sij,  As  before,  the  real  part  of  one  of  the 
n  —  r  columns  of  Q2  can  be  used  for 

MATRIX*  has  a  singular  value  decomposition  routine  called  SVD(matrix) 

[MAT86]. 


As  stated  previously,  once  c(J  has  been  selected,  it  is  important  to  check  the 
accuracy  of  the  approximation  of  Sij(s).  Since  A(s)  can  be  calculated  directly  at 
any  frequency  using,  for  example,  the  minimum  norm  solution,  it  is  important  to 
compare  those  values  against  the  calculation  of  using  the  coefficients  of  c,j. 


IV.  Selection  Criteria  for  the  QFT  Weighting  Matrix  A(s) 

Once  P(s)  and  P,(«)  are  known,  A(s)  can  be  calculated  using  the  Method 
of  Specified  Outputs  technique  of  Chapter  III.  Therefore,  the  key  to  determining 
A(s)  is  really  the  proper  choice  of  Pe.  An  equivalent  title  to  this  chapter  could  be 
Selection  Criteria  for  the  QFT  Equivalent  Plant  Matrix  P«(s). 

Selecting  an  acceptable  P*(j)  can  be  straightforward.  Each  element  of  the 
Pt{s)  matrix  is  a  transfer  function  that  relates  the  actual  system  outputs  Y  to 
the  desired  (or  commanded)  inputs  V  (F  is  also  the  equivalent  plant’s  input  com¬ 
mand  vector).  See  Figure  1.3.  The  transfer  function  is  a  “tracking”  function  or  a 
“rejecting”  function  depending  on  the  QFT  designer’s  wishes. 

Let  Pf  designate  a  failed  plant  matrix.  The  P/  matrix  is  determined  from 
Eq.  (2.6)  where  the  B  matrix  has  been  appropriately  scaled  in  order  to  model  a 
reduction  in  control  effectiveness.  Because  a  single  weighting  matrix  A(s)  is  used 
to  generate  the  equivalent  plant  matrix  Pt  throughout  the  entire  range  of  plant 
failures  Pf  ,  two  difficulties  may  arise.  First,  some  system  failures  are  catastrophic 
enough  that  the  failed  equivalent  plant  matrix  Ptf  =  P/A  no  longer  maintains  the 
necessary  or  desired  QFT  properties.  Second,  the  control  inputs  commanded  by 
A(s)  may  be  unreasonable  in  light  of  the  control  system’s  capabilities. 

The  following  discussion  contains  guidance  for  selecting  Pe{s)  which  in  turn 
determines  A(s)  when  using  the  Method  of  Specified  Outputs. 

J.l  Requirements  on  Pe(s) 

The  following  properties  of  Pe(s)  must  exist  if  QFT  is  to  be  used. 

4-1.1  Non-Singular  Pe(s)  For  Pe(s)  to  be  non-singular,  it  must  be  square 
and  its  determinant  non-zero.  That,  is,  P,(s)  must  be  controllable.  The  dimension 
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of  the  square  Pe  is  Z  x  /  where  l  is  the  number  of  outputs  that  the  designer  wishes  to 
control.  Note  that  l  <  m  where  m  is  the  number  of  independent  control  inputs. 

4-1.2  Equivalent  Plant  Element  Relationships  Depending  on  the  dimension 
of  the  equivalent  plant  matrix  Pe(s)  and  the  particular  QFT  MIMO  design  tech¬ 
nique  used,  a  specific  relationship  must  exist  between  the  elements  of  the  equivalent 
plant.  In  the  2x2  case,  the  diagonal  dominance  condition  must  hold.  Similar  con¬ 
ditions  exist  for  higher  dimensioned  equivalent  plant  matrices. 

4-1.3  Reasonable  Response  The  designer  should  choose  reasonable  transfer 
functions  for  the  elements  of  Pe(s).  If  the  output  is  to  follow  the  input,  then  a  first 
or  second  order  tracking  transfer  function  is  a  good  initial  choice. 

If  the  output  is  desired  to  be  zero  for  any  input,  then  a  first  or  second  order 
disturbance  rejection  transfer  function  is  a  good  initial  choice  for  the  elements  on 
the  main  diagonal  of  Pe(s)  [Lew88].  Either  disturbance  rejection  transfer  functions 
or  scalar  zeros  are  recommended  for  off-diagonal  elements  of  Pr(s).  See  the  cau¬ 
tionary  note  in  section  4.2.2  on  the  use  of  scalar  zeros  as  a  disturbance  rejection 
function  on  the  main  diagonal  of  Pe(s). 

To  determine  a  reasonable  choice  for  a  transfer  function  as  an  element  of 
Pf[s),  consider  the  control  input  hardware  limits.  Recall  that  this  analysis  assumes 
a  linear  system.  If  the  control  variables  becomes  saturated,  the  system  becomes 
non-linear. 

Another  consideration  in  selection  of  elements  of  Pe(s)  is  the  original  plant’s 
frequency  and  time  responses  to  the  standard  evaluation  inputs  in  each  of  the 
controls  ut.  For  example,  apply  a  step  input  of  Ui  to  the  original  plant  P  while 
holding  all  other  inputs  to  zero.  Observe  both  the  time  and  frequency  response  of 
system  outputs.  Do  not  ask  the  equivalent  plant  Pe(s)  to  provide  a  quantum  jump 
in  system  performance  when  compared  to  the  original  plant  P(s). 
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For  requirements  on  2x2  and  3x3  equivalent  plants,  see  section  2.4.1  of  Chap- 

1  ter  II. 

Computer  aided  design  tools  such  as  MATRIX*  can  only  simulate  the  re¬ 
sponse  of  proper  (i.e.  realizable)  transfer  functions.  Thus  the  degree  of  the  nu- 
^  merator  q  must  be  less  than  or  equal  to  the  degree  of  the  denominator  r.  The 

requirement  for  proper  transfer  functions  apply  not  only  for  the  elements  of  Pe{s) 
but  this  limitation  also  applies  to  elements  6,j(s)  of  the  weighting  matrix  A(s). 

►  Note  that  the  outputs  of  A (s)  are  the  control  vector  I[(s).  In  order  to  determine 
if  the  time  histories  u(t)’s  are  reasonable  for,  say,  step  inputs  to  A (s),  then  the 
^,j(s)’s  must  be  realizable. 

> 

4-2  Desired  Properties  of  Pe(s) 

The  following  properties  of  P,(s)  are  desired  in  order  that  the  QFT  design 
^  may  be  simplified. 

4-2.1  Minimum-Phase  <fcf[Pt(s)]  The  designer  should  select  the  zeros  of 
each  transfer  function  of  Pe{s)  so  that  det[Pe(s)]  are  all  minimum-phase  (i.e.  all 
t  zeros  of  the  determinant  lie  in  the  left-half  s-plane).  The  designer  must  check  to 

see  if  all  healthy  and  failed  equivalent  plant  determinants  are  minimum-phase  or 

if  the  right-half  s-plane  zeros  of  det[Pe/(s)j  are  so  far  into  the  r.h.p.  that  they  can 
^  be  ignored  [Hor88]. 

4-2.2  Diagonal  Pf{s)  A  diagonal  Pe(s)  is  desirable  so  that  non-diagonally 
related  inputs  and  outputs  are  decoupled.  A  note  of  caution  applies  when  using 
f  a  diagonal  Pe:  a  disturbance  rejection  transfer  function  on  the  main  diagonal  of  a 

Pe(s)  can  not  be  chosen  to  be  the  scalar  zero.  This  causes  det[Pf]  to  be  zero  and 

|  thus  violates  the  necessary  condition  1  of  section  2.4.1  and  the  MISO  Equivalent 

I 

)  method  can  not  be  applied  [Hou87:3-9j. 
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There  is  another  advantage  of  using  a  diagonal  Pt(s).  If  the  zeros  of  each 
transfer  function  of  a  diagonal  Pe(s)  are  chosen  to  be  minimum-phase,  then  the 
determinant  for  the  healthy  plant  will  necessarily  be  minimum-phase. 


V.  Example  of  the  Method  of  Specified  Outputs  Technique 


To  show  the  viability  of  the  Method  of  Specified  Outputs  technique  in  provid¬ 
ing  an  acceptable  equivalent  plant  for  QFT  design,  an  AFTI/F-16  aircraft  lateral- 
directional,  3-input /2-output  model  is  examined.  The  model  is  from  the  identical 
flight  condition  considered  in  [Sch86]. 

5.1  The  Model 

The  flight  condition  for  the  AFTI/F-16  is  Mach  0.9  at  20,000  ft.  The  block 
diagram  for  this  system  is  given  in  Figure  5.1.  See  Tables  5.1  and  5.2  for  system 
inputs  and  outputs. 


Ur  / 


«(/ 


Ur 


Figure  5.1.  AFTI/F-16  Lateral-Directional  3-Input  2-Output  Model 


Table  5.1.  Control  Inputs  and  Units 


Control  Inputs 

Description 

Units 

Ur  f 

right  flaperon  deflection 

degrees 

Ulf 

left  flaperon  deflection 

degrees 

Ur 

rudder  deflection 

degrees 

The  lateral-directional  state  space  model  is  given  in  Eqs.  (5.1)  and  (5.2) 


x  =  Ax  -f  Bu 


(5.1) 
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Table  5.2.  System  Outputs  and  Units 


Description 

Units 

sideslip  angle 

degrees 

roll  rate 

degrees  per  second 

y_  =  Cx  -f  Du  (5.2) 

where  x  is  the  vector  of  8  states,  u  is  the  vector  of  3  control  inputs,  and  t/  is  the 
vector  of  2  system  outputs.  See  Eqs.  (5.3)  thru  (5.5).  The  vector  of  state  variables 
is 


See  Table  5.3  for  a  description  of  the  8  states.  The  vector  of  control  inputs  is 

Table  5.3.  State  Variables  and  Units 


States  Description 


0  pitch  angle 


u  forward  velocity 


a  angle  of  attack 


pitch  rate 


bank  angle 


sideslip  angle 


Units 


degrees 


feet  per  second 


degrees 


degrees  per  second 


degrees 


roll  rate 

degrees  per  second 

yaw  rate 

degrees  per  second 

u  = 


(5.4) 


Urf 

Ulf 

Ur 


The  vector  of  system  outputs  is 


R  = 


0 

P 


The  state  space  model  plant  matrix  is 

A  =  [At  At  1 


where 


and 


0 

0 

0 

-32.183 

0.012075 

38.2906 

-0.00112 

-0.000022 

-1.48446 

-0.000309 

-0.00013 

4.27171 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0. 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0.03449 

-0.343554 

0.032636 

0 

-55.2526 

-2.80004 

0 

7.237 

-0.023184 

1 

-30.1376 

0.994789 

-0.777221 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-0.997556 

0.145674 

-0.36253 


(5.5) 


(5.6) 


(5.7) 


(5.8) 
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The  control  matrix  is 


0 

0 

0 

1.158405 

1.158405 

0 

-0.122462 

-0.122462 

0 

-3.236345 

-3.236345 

0 

0 

0 

0 

-0.0006855 

0.0006855 

0.037032 

-25.5251 

25.5251 

10.3955 

-0.62503 

0.62503 

-5.8089 

The  output  matrix  is 


C  = 


0  0  0  0  0  1  0  0 
0  0  0  0  0  0  1  0 


(5.9) 


(5.10) 


The  feedforward  matrix  is 

0  0  0 

(5.11) 

0  0  0 

The  AFTI/F-16  flaperon  deflection  limits  are  ±20  degrees  trailing  edge  up 
and  down  (TEU  and  TED).  The  rudder  deflection  limits  are  ±30  degrees  trailing 
edge  right  and  left  (TER  and  TEL).  The  flaperon  surface  rate  limit  is  52  degrees  per 
second  and  the  rudder  surface  rate  limit  is  120  degrees  per  second  [Bar83:Table  3]. 


5.2  The  Plant  Matrix,  Both  Healthy  and  Failed 

The  plant  matrix  relates  all  control  inputs  to  system  outputs.  For  the  exam¬ 
ple,  P(s)  has  the  dimensions  2x3,  i.e.  I  =  2  which  is  the  number  of  outputs  and 
m  =  3  which  is  the  number  of  inputs.  There  are  6  transfer  functions  contained  in 
P(s). 


P(s)  =  C'(sl-  A)~lB 


(5.12) 
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Define  the  healthy  plant  P  as  the  evaluation  of  Eq.  (5.12)  for  the  A,  B,  (7,  and 
D  matrices  of  Section  5.1.  Define  a  failed  plant  Pf  as  the  evaluation  of  Eq.  (5.12)  for 
the  A,C,  and  D  matrices  of  Section  5.1  and  for  a  scaled  B  matrix  that  incorporates 
a  reduction  in  control  effectiveness.  For  this  thesis,  the  healthy  plant  and  the  set 
of  failed  plants  of  Section  5.2.2  comprise  the  range  of  plant  uncertainty  that  this 
thesis  specifically  addresses. 


5.2.1  The  Healthy  Plant  For  the  8-state  model,  each  element  of  P  has  a 
common  denominator  polynomial  of  degree  8.  Upon  inspection  of  each  of  P(s)’s  6 
elements,  there  is  considerable  cancellation  between  the  numerator  zeros  and  the 
denominator  poles.  Therefore,  the  reduced  order  plant  matrix  is  represented  by 
6  numerator  polynomials  of  degree  3  and  one  denominator  polynomial  of  degree 
4.  The  time  histories  and  Bode  plots  of  the  original,  uncompensated  plant  are  in 
Figures  5.2  thru  5.7. 


(5.17) 

(5.18) 

(5.19) 


GAIN(DEG) 
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and 


dr  =  a*  +  3.5061a3  +  (1.1238e  +  l)aJ  +  (2.4367e  +  2 )s  +  (6.5450e  -  1)  (5.20) 

Note  the  symmetry  between  pn  and  pt2  and  between  p2i  and  p22.  The  pn  transfer 
function  describes  the  f3  response  to  a  right  flaperon  input  and  the  p\2  transfer 
function  describes  the  (3  response  to  a  left  flaperon  input.  For  a  healthy  aircraft, 
these  responses  are  equal  but  opposite  in  sign. 

The  poles  of  the  plant  matrix  (i.e.  the  roots  of  the  denominator  dr)  all  lie  in 
the  left-half  a-plane  and  thus  the  original,  uncompensated  plant  is  stable. 

-0.3910  + j'2.9616 
-0.3910  -  ;2.9616 

(5.21) 

-2.6969 
-0.02719 

5.2.2  The  Failed  Plants  As  stated  in  Section  5.2,  the  set  containing  the 
healthy  plant  and  all  the  failed  plants  comprise  the  range  of  plant  uncertainty 
for  this  thesis.  The  failures  considered  are  listed  in  Table  5.4.  An  assumption  is 
applied  where  plant  failures  can  be  modeled  by  reducing  the  effectiveness  of  the 
control  inputs  on  the  outputs  by  appropriately  scaling  the  B  matrix  of  Eq.  (5.12). 
For  example,  if  the  right  flaperon  suffers  a  failure  that  reduces  its  range  by  50 
percent,  then  the  healthy  B  matrix’s  first  column  is  multiplied  by  0.50.  Denoting 
this  failure  as  #  5,  then  Pja,iurt s  is  given  by  C(sl  -  A)~l  B/a,iure 5.  An  equivalent 
description  of  a  failure  is  illustrated  by  Eqs.  5.22  and  5.23  where  £',  is  the  positive 
scalar  less  than  or  equal  to  1  that  describes  the  effectiveness  of  the  ith  control 
surface. 

i  f  Urf 

E\Pn  E2P\2  E2P\2 

Vif 

E\P2\  e2p22  e3p23 

Ur 


(5.22) 
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Table  5.4.  Plant  Failures 


Failure  No. 

Surface  Effectiveness 

right  flaperon 

left  flaperon 

rudder 

1 

0% 

2 

0% 

3 

0% 

0% 

4 

0% 

0% 

5 

50% 

6 

50% 

7 

50% 

50% 

8 

50% 

50% 

9 

0% 

50% 

10 

50% 

0% 

11 

0.01% 

12 

0.01% 

13 

0.01% 

0.01% 

14 

0.01% 

0.01% 

15 

0.01% 

50% 

16 

50% 

0.01% 

17 

0.01% 

0.01% 

0.01% 

Surfaces  are  100%  Effective 

Jnless  Noted 

or 


Pn 

Pi  2 

Piz 

P2l 

P 22 

P2Z 

ExUrf 

E2Uif 

EiUr 


The  failures  in  Table  5.4  are  examined  for  this  example. 


(5.23) 


Note  that  the  stability  of  the  failed  plants  is  not  influenced  by  scaling  the 
B  matrix.  Since  P(s)  =  C’(sl  -  A)~l  B ,  then  only  the  common  denominator  due 
to  (si  —  _4)~l  determines  stability.  Since  this  thesis  assumes  that  the  stability 
derivatives  in  the  A  matrix  do  not  change  as  a  result  of  failures,  then  the  stability 
of  the  failed  plants  are  identical  to  the  healthy  plant. 

When  using  MATRIX*,  scaling  a  column  of  the  B  matrix  is  equivalent 


I 


to  scaling  the  appropriate  columns  of  NUM.  NUM  is  the  output  matrix  of  the 
[. NUM,  DEN ]  =  TFORM(S,NS)  routine.  When  using  MACSYMA,  scaling  a 
column  of  the  B  matrix  is  equivalent  to  scaling  a  column  of  the  plant  matrix  as  in 
Eq.  5.23.  See  Appendices  A  and  B. 

5.3  Selection  of  Equivalent  Plant 

An  equivalent  plant  matrix  Pe(a)  is  specified  by  the  QFT  designer  based 
on  the  considerations  of  Chapters  II  and  IV.  For  the  AFTI/F-16  example,  the 
equivalent  plant  is  a  2x2  matrix  of  transfer  functions  relating  the  2  inputs  (u!  and 
u2)  to  the  2  outputs  (sideslip  angle  (3  and  roll  rate  p). 

Recall  the  necessary  and  desired  characteristics  of  an  2x2  equivalent  plant 
for  MIMO  QFT  purposes.  The  matrix  Pt{»)  must  be  non-singular,  |peilpejj|  > 
|PenP«jil  as  u>  — ►  00  ,  det[Pt]  should  be  minimum  phase,  Pe  should  preferably  be 
diagonal,  and  the  output  of  P,  should  be  stable  and  reasonable  for  a  stable  P(s) 
and  reasonable  inputs. 

From  an  open-loop  performance  perspective,  the  goal  of  the  equivalent  plant 
is  to  reject  sideslip  and  follow  a  commanded  roll  rate.  Thus  peil  is  chosen  to  be 
a  disturbance  rejection  transfer  function,  pr„  is  chosen  to  be  a  tracking  transfer 
function,  and  both  pei3  and  pejl  are  chosen  to  be  the  scalar  zero. 

The  healthy  plant  P(s)  and  only  one  choice  for  Pe(s)  are  the  only  quantities 
needed  to  determine  a  weighing  matrix  A(s)  when  using  the  Method  of  Specified 
Outputs.  However,  in  order  to  compare  the  influence  of  the  choice  of  Pe(s)  on 
the  weighting  matrix  A(s),  several  Pe(s)'s  are  chosen.  See  Table  5.5.  A  common 
denominator  of  (5  4-  10)2  is  chosen  for  two  reasons.  First,  for  the  pej7{s)  roll 
rate  tracking  function,  this  corresponds  to  a  reasonable  roll  rate  response  for  an 
AFTI/F-16  type  aircraft.  Define  the  roll  mode  time  constant  Tr  as  63%  of  the 
time  it  takes  to  reach  a  steady  state  roll  rate  imparted  by  a  step  input  [Ros82:458]. 
With  this  definition,  Tr  is  0.22  seconds  as  determined  from  Figure  5.9  which  is 
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Table  5.5.  Choices  for  Pe(s) 


Choice 

Pe 

1 

(,+5)  n 

(.+10)»  U 

0  100 

U  (»+10)3 

2 

(<+9)  0 

(*+10)3  u 

0  100 

U  (e+10)3 

3 

(<+1)  0 

(t+10)2  u 

0  100 

u  (.+10)3 

4 

'  ('+50)  0 

('+10)*  u 

0  100 

U  (»+10)3  . 

5 

ton  0  ] 

(e+10)3  U 

0  100 

U  fe+10)3  . 

the  time  response  to  a  unit  step  input  to  the  transfer  function  Pe 22  =  . 

This  compares  favorably  with  Tn,„„r  of  1.0  seconds  for  a  class  IV  aircraft  during  a 
category  A  flight  phase  with  level  1  flying  qualities.  [Fly86:412j. 

Second,  in  order  to  compare  the  influence  of  the  numerators  of  peil  and  pr7J 
on  the  phase  of  the  determinant  of  the  failed  equivalent  plants  det[Pe/]  ,  then 
a  common  denominator  for  P(s)  is  a  reasonable  choice.  Note  that  for  all  the 
equivalent  plants  in  Table  5.5,  the  (3  rejection  transfer  function  peiI  does  not  have 
a  free  s  in  the  numerator  and  thus  a  steady  state  error  exists  for  step  inputs.  In 
light  of  the  desire  to  have  no  zeros  of  det[Pe]  or  det[Pe/]  in  the  right-half  s  plane  or 
on  the  imaginary  axis,  it  is  reasonable  to  achieve  f3  “rejection”  by  a  gain  much  less 
than  1  in  the  peil  transfer  function  and  not  by  a  free  s  in  the  numerator  of  ptu . 

Calculation  of  the  Weighting  Matrix  A (5) 

Recall  that  a  single  weighting  matrix  A(s)  is  used  to  transform  the  original 
plant  P{s)  into  the  equivalent  plant  Pe(s).  This  same  weighting  matrix  A (s) 
attempts  to  transform  a  failed  plant  P/(s)  into  the  equivalent  plant.  Of  course, 
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P/(a)A(a)  is  not  equal  P(a)A(a).  However,  it  is  hoped  that  Pe/(s)  =  Pf(s)A(s) 
will  maintain  all  the  necessary  and  desirable  QFT  characteristics. 

Since  P(a)  has  more  inputs  than  outputs  and  since  the  rows  of  P(s)  are  inde¬ 
pendent,  then  the  minimum  norm  solution  for  A(a)  applies.  Recall  the  minimum 
norm  solution  A(s)  =  P(s)T[P(s)P(s)T]_1Pe(s)  where  P(a)  is  the  healthy  plant. 
Since  P(s)  is  2x3  and  since  P*(s)  is  2x2,  then  A(s)  is  3x2. 

In  order  to  perform  a  simulation  of  the  response  of  A (a)  to  any  input,  each 
transfer  function  Sij(a)  of  A(s)  must  be  known.  The  outputs  of  A(s)  are  right 
flaperon  deflection  Urj(a ),  left  flaperon  deflection  Ui/(a),  and  rudder  deflection 
Ur(a).  Recall  that  the  coefficients  of  the  numerator  polynomial  and  denominator 
polynomial  for  6ij  are  represented  by  the  vector  c^.  Once  all  g^’s  are  found,  then 
a  simulation  of  the  frequency  and  time  response  of  each  Stj  is  conducted.  Each 
time  response  is  checked  for  reasonableness  in  light  of  the  AFTI/F-16  hardware 
limits. 

Each  of  the  6  transfer  functions  of  A(a)  are  calculated  using  the  numerical 
approach  to  the  Method  of  Specified  Outputs.  Since  P(s)  is  known  exactly  [i.e. 
no  measurements  are  needed  to  determine  the  frequency  response  of  P(j)],  then 
a  moderate  bandwidth  of  frequencies  is  used  to  calculate  A(s).  For  each  choice 
of  equivalent  plant,  A (ju>k)  is  numerically  evaluated  at  frequencies  of  w*  from  0.5 
to  15  radians  per  second  at  intervals  of  0.5  radians  per  second.  This  range  of  fre¬ 
quencies  is  in  the  neighborhood  of  the  bandwidth  of  the  original  plant’s  frequency 
response.  Note  that  /  =  30  is  the  number  of  frequencies  at  which  A(jw)  is  evalu¬ 
ated.  The  QFT  designer  chooses  the  maximum  degree  of  the  numerator  polynomial 
)  q  and  the  degree  of  the  denominator  polynomial  r  for  each  Sij.  The  30x[q  +  r  +  2j 

dimension  F  ntatrix  is  formed  for  6^  and  c,;  is  extracted  from  the  columns  of 
the  modified  Hermite  normal  form  matrix  F.  See  Appendix  A  for  an  example  of 
the  MATRIX.^'  command  files  used  to  numerically  determine  the  weighting  matrix 
transfer  functions 
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5.5  The  Failed  Equivalent  Plants 

For  this  example,  the  Method  of  Specified  Outputs  guarantees  that  the  prod¬ 
uct  of  the  healthy  plant  and  weighting  matrix  P(s)A(s)  is  exactly  the  specified 
healthy  equivalent  plant  Pe(s).  All  healthy  plants  are  selected  specifically  to  meet 
the  necessary  and  desired  QFT  characteristics. 

It  is  desired  that  all  equivalent  plants  in  the  range  of  plant  uncertainty  meet 
the  same  standards  as  the  healthy  equivalent  plant.  However,  if  a  failure  is  so  severe 
as  to  cause  the  failed  equivalent  plant  to  be  singular  and  thus  uncontrollable  or  to 
cause  det[Pe/]  to  be  non-minimum  phase,  then  the  QFT  designer  can  either  select 
another  P«(s)  and  generate  another  A(s)  or  the  designer  can  simply  exclude  that 
failure  from  the  set  of  failed  plants  and  proceed  with  the  QFT  design.  Also,  the 
QFT  designer  may  be  able  to  absorb  the  non-minimum  phase  case  in  the  design 
if  the  r.h.p.  zero  is  either  very  close  to  the  origin  or  reasonably  far  into  the  r.h.p. 
[Hor88]. 

For  this  example  problem,  the  determinant  of  the  failed  plants  det[Pe/]  is 
symbolically  evaluated  using  MACSYMA.  The  roots  of  the  numerator  polynomial 
of  det[P/A]  is  examined  for  right-half  s-plane  zeros  using  both  MACSYMA  and 
MATRIX*.  See  Appendix  B  for  an  example  of  the  MACSYMA  commands  used. 

5.6  Results 

5.6.1  The  A (s)  for  Each  Equivalent  Plant  A  frequency  sensitive  weighting 
matrix  is  determined  by  using  the  Method  of  Specified  Outputs.  For  each  choice  of 
equivalent  plant  in  Table  5.5,  a  weighting  matrix  A (s)  is  determined,  the  output  of 
A(s)  to  step  inputs  in  Vi(s)  and  V2(.s)  is  examined  for  reasonableness,  the  diagonal 
dominance  condition  is  checked,  and  det[Pe/j  is  examined  for  right-half  s-plane 
zeros. 
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L 


Pt  Choice  #  1 


For  Pe  = 


li±g] 

(*+10)3 

0 


0 

100 

(«+10)3 


,  the  weighting  matrix  is 


^nll 

&nl2 

6n2 1 

&n22 

^n31 

^r»32 

(6.0627e  -  5)a4  +  (1.1213e  -  2)s3  4-  (2.0606e  -  l)s2  +  a 

where 


(5.24) 


6„u  =  (3.3591  e  -  4)s4  4-  (1.8233e  -  3)a3  -  (7.1924e  -  3)j2  -  (3.9555e  -  2 )a  (5.25) 

Snl2  =  (  — 1.1966e  -  4)*3  -  (2.0199e  -  2)j2  -  (5.4795e  -  2 )a  -  (3.2570e  -  4)  (5.26) 

6nn  =  -*„n  (5.27) 

&n22  —  ~&n\2  (5.28) 

^n3i  =  (1.649Ge  -  3 )s*  +  (9.4161e  -  3 )aa  4-  (2.0183e  -  2).»2  4-  (7.1710e  -  2) a  (5.29) 
S„n  =  (— 4.4301e  -  6)s3  -  (1.3681e  -  3)s2  -  (1.9729e  -  4)a  -  (2.082 9e  -  3)  (5.30) 


Figures  5.8  and  5.9  show  the  desired  equivalent  plant  responses.  Fig¬ 
ures  5.10  and  5.11  show  the  time  histories  of  the  healthy  right  flaperon  and  healthy 
rudder  to  unit  step  inputs  of  1°  of  both  tq  and  t’j  into  the  weighting  matrix  A(.s). 
The  response  of  the  healthy  left  flaperon  is  identical  to  the  right  flaperon.  Since  the 
model  is  assumed  to  be  Unear,  then  total  right  flaperon  deflection  ur/  is  the  sum  of 
right  flaperon  deflection  due  to  a  1°  step  input  in  v\  and  right  flaperon  deflection 
due  to  a  1°  step  input  in  v2 .  The  statement  “the  total  response  is  the  sum  of  the 
individual  responses”  also  applies  to  the  rudder  deflection  response.  Recall  the 
hardware  Umits  on  the  AFTI/F-16  flaperon  and  rudder:  flaperon  deflection  is  20° 
TEU  and  TED;  rudder  deflection  is  30°  TER  and  TEL.  Based  upon  Figures  5.10 
and  5.11  and  the  control  surface  Umitations,  the  A(s)  in  Eq.  5.24  is  not  realistic. 
It  is  clear  that  inputs  greater  than  1°  saturate  the  control  ttr. 
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Figure  5.8.  Desired  Sideslip  Unit  Step  Response  for  Choice  #1 
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The  diagonal  dominance  condition  is  satisfied  for  all  failures  down  to  0.01% 
of  available  control  surface  deflections. 

Finally,  a  check  of  zeros  of  each  det[Pe/]  is  made.  Table  5.6  shows  the  zeros 
of  det(Pe/]  for  the  failures  of  Table  5.4.  In  summary,  this  choice  #1  of  Pe(s)  is  not 
acceptable  in  light  of  the  large  ur/,  u//,  and  uT  deflections  it  commands  for  unit 
step  inputs  of  vj  and  v 2. 


Pe  Choice  #2 
For  Pe  = 


('t9?  n 

(,+10)’  u 


0 


100 

(,+10)’  J 


,  the  weighting  matrix  is 


A  = 


where 


Snll  =  (3.3591e  -  4)34  +  (3.1669e  -  3)s3  -  (6.6175e  -  3)s*  -  (7.1199e  -  2)s  (5.32) 
Snl2  =  ( —  1.1966e  —  4)s3  —  (2.0199e  —  2)s2  -  (5.4795e  -  2)s  -  (3.2570e  -  4)  (5.33) 

6mi  =  ~6nn  (5.34) 

Sn22  =  -Sn  12  (5.35) 

^31  =  (1.1650e  -  3)/  +  (1.6401e  -  2)3  +  (2.4856e  -  2 )s2  +  (1.2908e  -  1)j  (5.36) 
Sn32  =  (— 4.4301e  -  6 )s3  -  (1.3682c  -  3)2  -  (1.9729c  -  5)3  -  (2.0829e  -  3)  (5.37) 


^nll  &nl2 

&n  21  ^n22 

^n31  ^n32 

(5  31 1 

(6.0627e  —  5)s4  4-  (1.1 

213e  -  2 )33 

+  (2.0606e-  l)32  +  3 

Figure  5.12  shows  desired  sideslip  f3  response  and  Figure  5.9  shows  the  desired 
roll  rate  p  response  for  unit  step  inputs.  Figure  5.13  is  similar  to  Figure  5.10  except 
that  the  right  flaperon  defection  has  slightly  less  overshoot  for  A(s)  choice  #2  than 
for  A(s)  choice  #1.  Figure  5.14  shows  less  rudder  deflection  than  Figure  5.11  but 
such  a  large  rudder  deflection  of  approximately  19°  for  1°  step  inputs  in  V'j  and  V2 
is  still  unacceptable. 


J 
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Table  5.6.  Zeros  of  det[Pe/]  for  Pe  Choice  #1 


Failure  No. 

Right-Half  Plane  Zeros 

Zeros  at  Origin 

1 

2 

2.9331 

0.0019  +  >0.0624 
0.0019  ->0.0624 

3 

2.7060 

0.4668  +  >4.0894 
0.4668  -  >4.0894 
0.0020  +  j 0.0624 
0.0020  -  j‘0.0624 

4 

5.3081 

1.9298  +>2.6871 
1.9298  -  >2.6871 

3 

5 

6 

7 

8 

9 

10 

4.3312 

0.6170 +  >4.6131 
0.6170  ->4.6131 
0.1036  +  >0.1609 
0.1036  ->0.1609 

11 

12 

13 

14 

15 

16 

17 

Table  5.7.  Zeros  of  det[Pe/]  for  P*  Choice  #2 


Failure  No. 

Right- Half  Plane  Zeros 

1 

2 

9.8535 

0.3633  +  >4.4728 
0.3633  -  >4.4728 

5 

3 

17.4698 

1.1118 

0.7381  +>5.1057 
0.7381  ->5.1057 

3 

4 

4.000 

1.4832 

3 

5 

6 

7 

8 

9 

10 

9.6381 

1.2592 

11 

12 

13 

14 

15 

16 

17 

The  diagonal  dominance  condition  is  satisfied  for  failures  down  to  0.01%  of 
available  control  surface  deflections. 

Table  5.7  shows  the  zeros  of  det[Pe/]  for  the  17  failure  cases. 

In  summary,  choice  #2  of  Pe(s)  is  not  acceptable  in  light  of  the  large  deflec¬ 
tions  of  the  control  surfaces  ur/,  Uj/,  and  ur.  Also,  the  “more”  minimum-phase 
det[P,(s)]  for  the  healthy  Pe(s)  actually  causes  the  zeros  of  det.[Pe/]  to  be  farther 
in  the  right-half  s  plane  than  choice  Pe(s)  #1. 
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’  ii±H 

For  P,(s)  =  (*+10,J 


0 


,  the  weighting  matrix  is 


L  0  u+w  J 


^nll 

f>nl2 

^n21 

&n22 

^n31 

^r»32 

(6.0627e  -  5)a4  +  (1.1213e  -  2)j3  +  (2.0606e  -  1  )s2  +  s 

where 


(5.38) 


6nll  =  (3.3591e  -  4)s4  +  (4.7963e  -  4 )a3  -  (7.7673e  -  3 )a2  -  (7.9110e  -  3)«  (5.39) 

Sn lt  =  (-1.1966*  -  4)s3  -  (2.0199e  -  2 )s2  -  (5.4795e  -  2)a  -  (3.2570e  -  4)  (5.40) 

Sn  si  =  -Sn  n  (5.41) 

^n22  =  —  f>nl2  (5.42) 

*„3i  =  (1.6496*  -  3)s4  +  (2.8178e  -  3)s3  +  (1.5510e  -  2)s2  +  (1.4342*  -  2)s  (5.43) 
K32  =  (-4.430U  -  6)s3  -  (1.3682*  -  2 )s2  ~  (1.9729e  -  4)s  -  (2.0829e  -  3)  (5.44) 

Figures  5.15  and  5.9  show  the  desired  equivalent  plant  responses  to  unit  step 
inputs.  The  right  flaperon  deflection  of  Figure  5.16  has  a  slightly  larger  overshoot 
than  Figures  5.10  and  5.13.  Figure  5.17  shows  an  excessively  large  rudder  deflection 
similar  to  Figure  5.11. 

The  diagonal  dominance  condition  is  satisfied  down  to  failures  of  0.01%  of 
available  surface  deflections. 

Table  5.8  shows  the  zeros  of  the  determinant  of  the  failed  equivalent  plant. 

In  summary,  choice  #3  for  Pt(s)  is  unacceptable  in  light  of  the  excessive 
control  deflections  commanded  by  A(s).  Also,  the  “less”  minimum-phase  healthy 
det[P,(.*)]  does  not  improve  the  location  of  the  r.h.p  zeros  of  det[Pf/(s)]  for  the 
failure  cases  involving  complete  surface  failures. 
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Table  5.8.  Zeros  of  det[Pe/j  for  Pe  Choice  #3 


Failure  No. 


Right- Half  Plane  Zeros 


Zeros  at  Origin 


1 


8 

0.3956 

0.3956 


0687 

+  >1.2675 
->  1.22675 


6, 

0.3657 

0.3657 


6973 

+  >1.5310 
->1.5310 


5 

0.9213 

0.9213 

0.3507 


,2977 
+  >3.1922 
-  V3.1922 
+  '0.9407 


5 

6 

7 

8 
9 


10 


11 

12 

13 

14 

15 

16 
17 


0.3507  -  >0.9407 


30.9174 

1.2606 

0.4766 

4.4046  +  >28.9746 
4.4046  -  >28.9746 
0.1071  +>0.4238 
0.1071  ->0.4238 


For  Pt  Choice  #4 
For  Pe  = 


(«+S0) 

(«+10)J 


0 


0 

100 
(«  +  10)J  J 


,  the  weighting  matrix  is 


A  = 


^nll  i>n  12 

&n2l  ^r»22 

^n31  Sn32 

E _  - 

(6.0627c  —  5).»4  +  (1.1 

1213c  -  2)s3  -1-  (2.0606e  -  l)s2  +  s 

(5.45) 


where 


6nll  =  (3.3591e  -  4 )s*  +  (1.6939e  -  2 )s3  -  (7.24910e  -  4 )s2  -  (3.9555e  -  l)a  (5.46) 
Sni2  =  (  — 1.1966e  -  4)s3  -  (2.0199e  -  2)s2  -  (5.4795e  -  2 )s  -  (3.2570e  -  4)  (5.47) 

^n21  =  —  ^nll  (5.48) 

^n22  =  —  £t»12  (5.49) 

Sn3l  =  (1.6496e  -  3)s<  +  (8.3648e  -  2)j3  +  (7.2751e  -  2)s2  +  (7.1710e  -  l)a  (5.50) 
*n32  =  (~4.4301e  -  6)s3  -  (1.3682e  -  3 )a*  ~  (1.9729c  -  4 )s  -  (2.0829e  -  3)  (5.51) 


Figures  5.18  and  5.9  show  the  desired  equivalent  plant  responses  for  unit  step 
inputs  in  iq  and  vj-  The  right  flaperon  deflection  of  Figure  5.19  does  not  reach  0° 
as  quickly  as  Figures  5.10,  5.13,  and  5.16.  The  rudder  deflection  of  Figure  5.20  is 
similar  to  that  of  Figure  5.11  except  that  ur  of  Figure  5.20  does  not  reach  0°  as 
quickly. 

The  diagonal  dominance  condition  is  satisfied  for  all  failures  other  than  ones 
involving  100%  surface  failure. 

Table  5.9  shows  the  zeros  of  det[Pe/]  for  the  failures  considered. 

In  summary,  Pe(s)  choice  #4  is  unacceptable  in  light  ofjl.e  large  surface 
deflections  commanded  by  A(s).  Also,  the  presence  of  (s-f-50)  in  ti  e  numerator  of 
p,u(s)  does  not  appreciably  improve  the  location  of  the  r.h.p.  zero*  of  det[P^(.s)] 
for  the  failures  considered. 
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RUDDER  (DEG) 


Figure  5.20.  Rudder  Unit  Step  Response  for  Choice  #4 


Table  5.9.  Zeros  of  det[Pe/j  for  Pt  Choice  #4 


For  Pt  Choice  #5 


For  Pt  = 


0.01 

(*+10)3 

0 


0 

100 

(»+10)3 


,  the  weighting  matrix  is 


^nll  &nl2 

&n21  &n22 

^n31  f>n3  2 

(6.0627e-5)s4  +  (l.l 

213e-2)s3 

+  (2.0606e  -  l)s2  +  3 

where 

Snn  =  (3.3591e  -  6 )s3  +  (1.4372e  -  6 )s2  -  (7.9110e  -  5 )s  (5.53) 

*n!2  =  (-1.1966C  -  4)*3  -  (2.0199e  -  2 )s2  -  (5.4795e  -  2 )s  -  (3.2570e  -  4)  (5.54) 

*„21  =  Snu  (5.55) 

*n22  =  -6„n  (5.56) 

6„31  =  (1.6496e  -  5)s3  +  (1.1682e  -  5)s2  +  (1.4342e  -  4)a  (5.57) 

Sn32  =  (-4.430U  -  6)s3  -  (1.3682e  -  3)a2  -  (1.9729e  -  4 )s  -  (2.0829)  (5.58) 


Figures  5.21  and  5.9  show  the  desired  equivalent  plant  responses  to  unit  step 
inputs  in  t’j  and  V2-  Figures  5.22  and  5.23  show  the  right  flaperon  deflection  for 
unit  step  inputs  in  tq  and  iq.  Note  the  reasonable  magnitude  of  urf  for  the  1°  step 
inputs.  Figures  5.24  and  5.25  show  the  rudder  deflection  for  unit  step  inputs  in  tq 
and  iq.  Note  the  reasonable  magnitude  of  ur  for  the  1°  step  inputs.  In  terms  of 
reasonable  response  for  reasonable  inputs,  choice  #5  for  Pt(s )  is  acceptable. 

All  failed  equivalent  plants  are  non-singular  for  all  failures  not  involving  0% 
surface  failures.  Also,  the  diagonal  dominance  condition  is  satisfied  for  all  failures 
not  involving  0%  surface  failures. 

Table  5.10  shows  the  zeros  of  det[Pt/(s)].  Note  that  det[Pe(s)]  for  the  healthy 
Pe{s)  has  no  zeros  whatsoever.  Yet,  the  failures  involving  100%  surface  failures 
continue  to  have  r.h.p  zeros. 
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5 

TIME(SEC) 


6 


8 


Desired  Sideslip  Unit  Step  Response  for  Choice  #5 
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RIGHT  FLAPERON  (DEG) 
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Table  5.10.  Zeros  of  det[/,e/(j)j  for  Pe(.»)  Choice  #5 
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Based  on  the  above,  the  weighting  matrix  of  Eq.  (5.52)  produces  acceptable 
results  for  failures  up  to  99.99%  of  available  control  surface  deflections.  For  com¬ 
parison  with  the  healthy  equivalent  plant  unit  step  response  Figures  5.26  and  5.27, 
the  unit  step  responses  of  failed  equivalent  plants  #12,  #13,  and  #14  are  shown 
in  Figures  5.28  thru  5.33.  Recall  the  failure  cases  are  listed  in  Table  5.4. 

5.6.2  Summary  All  the  choices  for  Pe{s)  provide  some  open-loop  perfor¬ 
mance  improvement  over  the  uncompensated  plant  P{s)  of  the  lateral-directional 
model  of  the  AFTI/F-16  at  Mach  0.9  at  20,000  feet.  For  all  healthy  plants,  the 
weighting  matrix  generated  by  the  Method  of  Specified  Outputs  provides  the  nec¬ 
essary  and  desired  structure  to  the  equivalent  plant  for  the  QFT  design  process 
to  proceed.  The  desired  characteristic  of  minimum  phase  det[Pe/(s)]  for  all  plant 
failures  can  not  be  met  with  the  choices  of  Pe(j)  considered.  However,  failures 
up  to  0.01%  of  available  control  surface  deflection  are  tolerable  when  using  the 
weighting  matrix  A(s)  generated  from  the  equivalent  plant  matrix  Pr(.s)  of  choice 
#5  of  Table  5.5.  Thus,  for  failures  up  to  0.01%  of  available  control  authority,  the 
Method  of  Specified  Outputs  generates  a  frequency  sensitive  weighting  matrix  A(a) 
that  transformed  the  original  plant  matrix  P(s)  into  and  equivalent  plant  matrix 
Pt(s)  that  possesses  the  necessary  and  desired  characteristics  for  the  QFT  design 
to  proceed. 

Note  that  failures  2,3,4,  and  10  of  Table  5.4  cause  the  det [Pe  ]’s  to  be  non¬ 
minimum  phase  for  all  5  choices  of  Pt(s)  .  Also,  these  failures  cause  Pe/(s)  to  be 
almost  singular.  Failures  2,3  and  10  involve  complete  loss  of  the  rudder.  Failure 
4  involves  complete  loss  of  both  flaperons.  On  the  other  hand  failure  1,  which 
involves  the  loss  of  only  one  flaperon  while  the  other  independent  ftaperon  remains 
healthy,  results  in  all  det[Pe/]’s  having  no  r.h.p.  zeros  and  is  thus  minimum-phase 
for  all  5  choices  of  Pe(s).  Further  exploration  is  required  to  find  an  equivalent  plant 
that  will  make  all  det[Pe/]’s  minimum-phase. 


5-46 


ACTUAL  P  (DEG/S) 


FAILURE  12  BETA  (DEG) 


FAILURE  12  P  (DEG/S) 


FAILURE  13  BETA  (DEG) 
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Figure  5.31.  Roll  Rate  Unit  Step  Response  for  Pe/  Failure  #13 

i 
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FAILURE  14  P  (DEG/S) 


VI. 


Conclusions  and  Recommendations 


6.1  Conclusions 

The  Method  of  Specified  Outputs  is  successful  in  generating  a  frequency 
sensitive  weighting  matrix  A (s)  for  use  in  the  MIMO  design  process.  However,  in¬ 
suring  that  A(.j)  commands  reasonable  control  variables,  insuring  that  the  diagonal 
dominance  condition  is  met,  and  insuring  that  all  det[P«/(s)]  are  minimum  phase 
still  requires  a  trial-and-error  approach.  The  advantage  of  this  trial- and-error  ap¬ 
proach,  when  compared  to  earlier  weighting  matrix  selection  techniques,  is  that 
choices  for  the  equivalent  plant  constitutes  the  trial-and-error  and  not  the  choices 
of  the  weighting  matrix  itself.  In  cases  where  the  number  of  system  outputs  is  less 
than  or  equal  to  the  number  of  control  inputs,  choosing  a  desired  equivalent  plant 
with  the  necessary  and  desired  QFT  characteristics  does  guarantee  the  actual, 
healthy  equivalent  plant  will  possess  the  necessary  and  desired  Q FT  properties. 

For  the  AFTI/F-16  example  considered,  a  frequency  sensitive  weighting  ma¬ 
trix  A(s)  is  found  that  transforms  the  original  plant  matrix  P(s)  into  an  equivalent 
plant  matrix  P«(s)  that  satisfies  all  QFT  necessary  and  desired  characteristics  for 
all  failures  up  to  a  99.99%  loss  in  control  surface  deflection  capability.  Complete 
restriction  of  travel  of  either  the  rudder  or  both  flaperons  results  in  failed  equiva¬ 
lent  plants  that  do  not  satisfy  the  QFT  requirements.  Once  the  transfer  function 
of  the  weighting  matrix  has  been  determined  numerically,  use  of  a  symbolic  math¬ 
ematical  tool  like  MACSYMA  allows  the  QFT  designer  to  quickly  compare  the 
affects  of  different  weighting  matrices  on  the  original  plant  matrix. 

6.2  Recommendations 

Future  investigation  in  this  area  should  include: 
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•  Investigate  the  use  of  MACSYMA  for  determining  A (s)  from  the  minimum 
norm,  least  squares,  and  exact  square  solutions. 

•  Consider  the  determination  of  A (s)  as  an  optimal  control  problem  with  path 
constraints  and  not  just  end-point  constraints.  Thus  the  controls  commanded 
by  A(s)  should  satisfy  the  hardware  limitations  of  the  control  system  for 
reasonable  desired  system  outputs. 

•  A  systems  identification  scheme  such  as  LTC  Baker’s  should  provide  the 
determination  of  the  coefficients  of  A (s)  even  with  the  presents  of  noise  in 
the  measurement  of  the  plant  matrix  P(s).  Also,  LTC  Baker’s  technique 
explicitly  provides  a  common  denominator  for  the  weighting  matrix  A(s). 
[Bak80], 

•  Consider  the  role  of  the  weighting  matrix  as  it  steers  the  eigenvectors  of  P 
through  Pt  space.  The  frequency  sensitive  weighting  matrix  provides  fre¬ 
quency  weighted  steering  through  the  input  space  to  the  output  space.  Thus 
A(s)  gives  the  QFT  designer  control  over  which  eigenvectors  are  excited  and 
at  which  frequency  they  are  excited  [Lew88]. 

•  If  the  QFT  designer  wishes  to  use  a  weighting  matrix  of  constant  elements 
instead  of  frequency  sensitive  elements,  then  Minimum  Control  Effort  anal¬ 
ysis  might  be  used  to  determine  weightings  for  the  constant  elements  6tj ’s 
[Rei83:387]  and  [May87]. 


6-2 


Appendix  A.  MATRIXx  Command  Files 


Examples  of  MATRIX*  command  files  are  presented  in  this  appendix  to 
illustrate  how  to  input  a  state  space  model,  determine  the  transfer  functions  that 
reside  in  the  original  plant  matrix  P(s)  and  the  failed  plants  P/(s),  and  use  the 
Method  of  Specified  Outputs  to  numerically  determine  the  weighting  matrix  A(s) 
that  transforms  P(s)  into  the  desired  equivalent  plant  Pe(s)-  Only  skeletons  of 
vital  sections  of  the  working  command  files  are  included. 

//MATRIXx  Command  File 

/ /This  program  is  intended  to  determine  the  original  plant’s  transfer 
//functions,  both  healthy  and  failed,  and  to  determine  if  the 
//healthy  and  failed  plants  are  controllable 
A  =  [state  space  model  plant  matrix] 

B  =  [state  space  model  control  matrix] 

C=(state  space  model  output  matrix] 

D  =  [state  space  model  feedforward  matrix] 
ns = number  of  states 
l=number  of  outputs 
m= number  of  inputs 


A-l 


program  continued 

//head thy  plant  matrix  of  transfer  functions 
smat=[A  B;C  D] ; 

[num,den]  =tform(smat,ns) 

/ /plant  matrix  transfer  functions  with  a  50%  failure 
/ /of  the  first  control  input 
Bfail=B;Bfail(:,l)=B(:,l)*.5;sBfail=[A  BfailjC  D]; 
[numBfail,den]  =tform(sBfail,ns) 

//ask:  is  the  healthy  system  controllable? 
[scon,nscon]  =cntrlable(smat,ns) 

//ask:  is  the  failed  system  controllable? 

[scon,nsconf]  =cntrlable(sBfail,ns) 


//MATRIX*  Command  File 

/ /This  program  is  intended  to  determine  the  numerical  value 
//of  the  weighting  matrix  DELTA(jw)  evaluated  at  specific  frequencies. 
//DELTAMAT  is  the  matrix  whose  first  m  rows  is  the  DELTA(jw)  matrix 
//evaluated  at  the  first  frequency,  its  second  m  rows  is  the  DELTA(jw) 
//matrix  evaluated  at  the  second  frequency,  ect. 

/ /f=maximum  frequency  at  which  DELTA(jw)  is  evaluated 
//m= number  of  control  inputs;  for  this  example,  let  m=3 
//l=number  of  system  outputs;  for  this  example,  let  1=2 
/ /ns=number  of  states;  for  this  example,  let  ns=4 
//w=frequency  in  rad/s 

//pd=denominator  of  DELTA(jw)  matrix  evaluated  at  frequency  w 
//pn(i  j)=numerator  of  element  ij  of  DELTA(jw)  matrix  evaluated  at  w 
//patw=evaluation  of  the  lxm  plant  matrix  at  frequency  w 
//pmat=matrix  composed  of  patw  matrices  evaluated  at  different  w’s 
//Pemat=desired  plant  matrix  specified  by  the  QFT  designer 
//for  this  example,  let  Pemat  =  [(.01)/((s+10)**2),0;0,100/((s+10)**2)] 
//den=matrix  from  previous  example 
//num=matrix  from  previous  example 


I 


1 


A -3 


program  continued 
for  w=l:f,... 
s=jay*w;... 
pd=0; 

for  i=l:(ns+l),... 
pd=pd+den(i)*s**(ns-(i-l ));... 
end;... 

pn(l,l)=num(l,4)*s**3+num(l,7)*s**2+num(l,10)*s+num(l,13);... 

pn(l,2)=num(l,5)*s**3+num(l,8)*s**2+num(l,ll)*s+num(l,14);... 

pn(l,3)=num(l,6)*s**3+num(l,9)*s**2+num(l,12)*s+num(l,15);... 

pn(2,l)=num(2,4)*s**3+num(2,7)*s**2+num(2,10)*s+num(2113);... 

pn(2,2)=num(2,5)*s**3+num(2,8)V*2+num(2,ll)*s-i-num(2,14);... 

pn(2^)=num(2,6)*s**3-fnum(2,9)*s**2+num(2,12)*s+num(2,15);... 

for  i=l:l,... 

for  j=l:m,... 

patw(i,j)=pn(ij)/pd;... 

end;... 

end;... 

pmat([(w-l)*l+l:w*l],:)=patw;... 

Pemat([(w-l)*l+l:w*I],:)=[(.01)/((s+10)#1*2),  0;  0,  100/((s+10)**2)] 
DELTAMAT([(w-l)*m+l:w*m], :)=... 

patw’*inv(patw*patw’)*Pemat(((w-l)*l-l-l:w*l] 

DELTAatW=DELTAMAT([(w-l)*m+l:w*m], :);... 


//MATRIXx  Command  File 

//This  program  is  intended  to  determine  the  coefficients  of 
/ /element  ij  of  the  weighting  matrix  DELTA(jw). 

//This  program  does  not  place  the  -1  in  the  appropriate 
/  / column  of  the  row  reduced  echelon  matrix;  the  user  must 
/ /do  that 

//deltai=row  number  of  DELTA(jw)  element 
//deltaj=column  number  of  DELTA(jw)  element 
// q=degree  of  numerator  polynomial  of  DELTA(jw)  element 
/ /r= degree  of  denominator  polynomial  of  DELTA(jw)  element 


A-5 


program  continued 
for  w=l:f,... 
s=jay*w,... 

eldelij(w)=DELTAMAT((w-l)*m+deltai,deltaj);... 

LittleN(w,:)=-[s**q,  s**(q-l),  s**(q-2),  s**(q-3),  s**(q-4),  s**(q-5),... 
s**(q-6),  s**(q-7),  s**(q-8),  s**(q-9),  s**(q-10)];... 
BigN(w,:)=LittleN(w,[l:q+l]);... 

TinyD(w,:)=(s**r,  s**(r-l),  s**(r-2),  s**(r-3),  s**(r-4),  s**(r-5),... 
s**(r-6),  s**(r-7),  s**(r-8),  s**(r-9),  «**(r-10)];... 
LittleD(w,:)=TinyD(w,[l:r+l]);-.. 

BigD(w,:)=eldelij(w)*LittleD(w, 

end; 

BigF =[BigN  ,-BigD] 

sizeBigF =size(  BigF  ),rankBigF =rank(  BigF ) 
ech=rref(BigF) 

//The  goal  of  this  program  is  to  find  any  vector  of  real 
/ /coefficients  that  lie  in  the  nullspace  of  the  BigF  matrix. 

/ /If  BigF  is  not  singular  then  no  nullspace  exists  and  the  user 
//must  select  higher  values  of  q  and  (or)  r  so  that  BigF  is  rank 
//deficient 

//Again,  the  user  must  replace  the  main  diagonal  0  with  a  -1  from  any 
/ /column  of  the  row  reduced  echelon  BigF  matrix  and  take  the  real 
//part  of  that  column  vector.  The  first  [q+r+2]  elements  of  that 
//modified  column  of  the  row  reduced  echelon  BigF  matrix  (cij)  lies 
/ /in  the  nullspace  of  BigF.  As  a  check,  evaluate  the  product  BigF*cij 
//to  see  if  the  elements  of  the  resulting  column  vector  is  almost 
//the  zero  vector  (e.g  10E-10  magnitude  for  each  element). 

//Once  all  the  transfer  functions  of  the  weighting  matrix  are 
//determined,  it  is  advantageous  to  use  MAC'SYMA  to  perform  the 
//matrix  multiplication  that  transforms  P(s)  into  Pe(s). 
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The  following  is  a  summary  of  MATRIX*  built-in  functions  that  the  user 
may  find  useful.  They  include  CNTRLABLE,  FOR  LOOPS,  INV,  JAY,  POLY, 
RANK,  REAL,  ROOTS,  RREF,  SIZE,  SVD,  and  TFORM.  , 
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Appendix  B.  Useful  MACSYMA  Commands 


This  appendix  provides  an  example  of  MACSYMA  commands  for  manipu¬ 
lation  of  the  matrix  product  Pe(s)  =  P(s)A(s).  Let  the  original  system  have  3 
inputs  and  2  outputs  and  let  the  original  plant  matrix  be  as  given  in  Eq.(B.l)  and 
let  let  the  weighting  matrix  be  as  given  in  Eq.  (B.2). 

(a  +  1)  (a  4-  2)  (a  4-  3) 

(a  4  4)  (a  +  5)  (a  +  6) 


P(a)  = 


A(a)  = 


(a  +  10) 

(a  4-  0.1)  (a  4  0.2) 
(a +  0.3)  (a  +  0.4) 
(a  +  0.5)  (a  +  0.6) 


(B.l) 


(B.2) 


(a  +  0.1) 

The  following  MACSYMA  commands  show  how  to  input  the  P(a)  and  A(a) 
matrices,  determine  the  actual  healthy  and  failed  equivalent  plant  matrices  for  a 
50%  loss  in  control  input  #1,  find  the  determinant  of  the  actual  equivalent  plants, 
and  find  the  zeros  of  the  determinants.  In  the  following,  Pn  is  the  matrix  of  nu¬ 
merator  polynomials  of  P(a),  Pd  is  the  common  denominator  polynomial  of  P(a), 
Pnf  is  the  matrix  of  numerator  polynomials  of  P(s)  where  control  input  #1  has 
lost  50%  of  available  control  effectiveness,  Dn  is  the  matrix  of  numerator  polyno¬ 
mial  of  A(s),  and  Dd  is  the  common  denominator  polynomial  of  A(s).  Also,  Pen 
is  the  matrix  of  healthy  equivalent  plant  numerator  polynomials,  Ped  is  the  com¬ 
mon  denominator  polynomial  of  Pe(s),  and  Penf  is  the  matrix  of  failed  equivalent 
plant  numerator  polynomials.  Finally,  detPen  is  the  determinant  of  the  healthy 
equivalent  plant,  detPenf  is  the  determinant  of  the  failed  equivalent  plant,  zdet- 
Pen  are  the  zeros  of  the  healthy  equivalent  plant  determinant,  and  zdetPenf  are 
the  zeros  of  the  failed  equivalent  plant  determinant.  Recall  that  Pe(s)  -  P(s)A(s). 
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Pn:matrix([(«  v  l),(s+2),(s+3)],[(s+4),(s+5),(s+6)]); 
Pd:(s+10); 

Puf(l,l]:0.5*Pn[l,l]; 

Pm{2,l]:0.5*Pn[2,l]; 

Dn:matrix([(s-f0.1),(s+0.2),(s+0.3)J1[(s+0.4),(s+0.5),(s4-0.6)]); 

Dd:(s+0.1); 

Pen-.expand(Pn.Dn); 

Ped:expand(Pd*Dd); 

Penf:expand(Pnf.Dn); 

detPen:expand(determinant(Pen)); 

detPen:expand(Pen[l,l]*Pen[2,2]-Pen[l,2]*Pen[2,lj); 

det  Penf:expand(  determinant  ( Penf ) ) ; 

zdetPen:al!roots(detPen); 

zdetPenf:ailroots(detPenf); 


Other  useful  MACSYMA  commands  are  writefile( filename),  closefile(), 
save( filename),  load(filename),  %  ,  and  quit(). 
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Ada88. 

Bak80. 

Bar83. 

Fly  86. 
Gel74. 
Ham87. 

IIor79. 

Hor88. 

Hou87. 

Lan69. 

Lev59. 

Lew88. 

MAT86 

May87. 
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